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Abstract. We use the worldline representation of field theory together with 
a variational approximation to determine the lowest bound state in the scalar 
Wick-Cutkosky model where two equal-mass constituents interact via the ex- 
change of mesons. Self-energy and vertex corrections are included approx- 
imately in a consistent way as well as crossed diagrams. Only vacuum- 
polarization effects of the heavy particles are neglected. In a path integral 
description of an appropriate current-current correlator an effective, retarded 
action is obtained by integrating out the meson field. As in the polaron prob- 
lem we employ a quadratic trial action with variational functions to describe 
retardation and binding effects through multiple meson exchange. The vari- 
ational equations for these functions are derived, discussed qualitatively and 
solved numerically. We compare our results with the ones from traditional ap- 
proaches based on the Bethe-Salpeter equation and find an enhanced binding 
contrary to some claims in the literature. For weak coupling this is worked out 
analytically and compared with results from effective field theories. However, 
the well-known instability of the model, which usually is ignored, now ap- 
pears at smaller coupling constants than in the one-body case and even when 
self-energy and vertex corrections are turned off. This induced instability is 
investigated analytically and the width of the bound state above the critical 
coupling is estimated. 



1 Introduction 

The bound state problem in Quantum Field Theory has a long history start- 
ing with the classic paper by Salpeter and Bethe more than 50 years ago. 
While the Bethe-Salpeter equation (BSe) is formally exact, its building blocks 
are only defined as infinite series of increasingly higher n-point Green functions. 
This hierarchy of coupled equations (called Dyson-Schwinger equations) has to 
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be truncated for practical purposes. In this way approximations are introduced 
which are hard to control and sometimes violate symmetries (e.g. gauge invari- 
ance). The most common approximation is the ladder approximation where bare 
propagators and vertices are used. There is a vast literature on the ladder BSe 
and its many deficiencies (e.g. existence of abnormal states, failure in the large 
mass limit etc.) Practical methods to include an extended set of diagrams 
(dressed propagators 3 or cross-ladder [3]) have been developed only recently. 

Many other variants have been investigated over the years: 3-dimensional re- 
ductions (also called the "quasipotential" approach) [SIIEIIZI, relativistic Hamil- 
tonian models in light-front [S] and point form 9 quantization, Feynman- 
Schwinger representations JO] to name only a few. It is close to impossible (and 
not intended here) to give a comprehensive overview over all attempts and the- 
oretical approaches available in the literature. 

Apart from the theoretical interest, there are (at least) three areas of physics 
in which the relativistic bound state problem needs to be addressed: the first 
one is bound-state Quantum Electrodynamics (QED) where ultra-precise exper- 
imental data are available. However, from the theoretical point of view this is 
a special field since one can employ perturbation theory in powers of the fine- 
structure constant and its logarithm after binding between the constituents has 
been established. Impressive accuracy matching the experimental precision has 
been achieved in recent years jllj . 

There is another area for relativistic bound-state calculations where strong- 
coupling is required: the electromagnetic structure of nuclear few-body systems 
at intermediate and high energies. For example, relativistic effects in the elec- 
tromagnetic form factor of the deuteron seem to be important already at sur- 
prisingly low momentum transfers j!21 ll.'jj . Unfortunately, these effects cannot 
be clearly separated from those of the phenomenological nuclear interaction and 
of the assumed form of the current operator. 

Finally there is the area of hadronic physics where one has to understand 
how (light) relativistic quarks and gluons bind to form the low-energy mesons 
and baryons (possibly also hybrids and exotics). Here, Quantum Chromodynam- 
ics (QCD) at least provides a clear underlying theoretical framework. However, 
the problems of confinement, chiral symmetry breaking and strong coupling are 
formidable and have not been solved (analytically) at low energies. Lattice Gauge 
Theory (see, e.g. ref. |14j ) is considered as the prime method to obtain gauge- 
invariant results from first principles, albeit with enormous numerical effort and 
problems of its own. The continous progress of lattice calculations not withstand- 
ing, considerable progress has also been made in the last years in solving the 
Dyson- Schwinger equations for Landau-gauge QCD under some simplifications 
|15j and in describing the low-lying hadrons as bound-state of quarks and gluons 
[T5| . While phenomenologically quite succesful and often going beyond the ladder 
BSe, these calculations still have limitations due to truncations, gauge depen- 
dence and the use of model propagators. In a similar way as in QED, remarkable 
theoretical progress has been achieved recently for special systems like heavy 
quarkonia by using the techniques of (nonrelativistic) effective field theories ^Fj\ . 

Given the twisted history of the relativistic bound-state problem, it might be 
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useful to consider it from a different angle, namely from the one provided by the 
worldline variational approach. Variational methods are widely used in molecu- 
lar, atomic and nuclear physics as they offer accurate, non-perturbative results 
for bound and scattering states in a quantum-mechanical framework. In field the- 
ory their use is more restricted |18j since one cannot handle non-gaussian wave 
functionals as trial states. Nevertheless, there have been a number of variational 
calculations for relativistic bound-state systems, most notably by Darewych and 
collaborators |191 I2U| . In contrast to these works we advocate here the world- 
line description of field theory which entails a huge reduction of the dynamical 
degrees of freedom. This is essential for a variational calculation in which only 
quadratic trial actions can be employed for analytical evaluations. The clas- 
sic example is Feynman's treatment of the polaron ~~ a nonrelativistic field 
theory for an electron interacting with the phonons in a crystal - which still 
stands out among all approaches describing this system non-perturbatively for 
all coupling constants. The key to this success was that, in a path integral treat- 
ment of the problem, he could integrate out the phonons exactly leading to a 
retarded two-time action for the electron. For the latter a variational approxi- 
mation then gave results unmatched by all other (analytical) methods. Shortly 
after this seminal paper, Mano, a student of Feynman's, applied these results 
to a scalar relativistic theory [22]. However, his pioneering work was largely 
forgotten as was the description of relativistic systems by trajectories (world- 
lines) . Recent times have seen a revival of worldline methods [22] under different 
names (e.g. Feynman-Schwinger representation) and for different purposes, e.g. 
for an efficient calculation of diagrams or effective actions. In a series of papers 
(OESlEniEZlEHlEn]) we have extended the "polaron" variational approach 
to field theory as pioneered by Mano. This was first done in the context of the 
scalar Wick-Cutkosky (WC) model [20] and later extended to a more realistic 
fermionic theory, viz. QED f |311 I32| ). As an attempt to catch non-perturbative 
features of these field theories, the results have been encouraging: for example, 
the well-known instability of the WC model [33] was detected by the worldline 
variational approach, and in QED a fully gauge-invariant approximate result for 
the anomalous mass dimension of the electron was obtained |34j . 

Therefore we believe that a fresh look at the relativistic bound-state problem 
from the worldline variational perspective may be useful even in an unrealistic 
model like the WC model, which incidentally was developed just for the purpose 
of studying bound states . In previous applications [21] - [2H] we used a version 
of that scalar model which describes neutral "nucleon" fields interacting with 

"meson" fields x- It is given by the Lagrangian (our metric signature is H 

and in general we put H = c = 1). 

£>l = \ " \M 2 ^ 2 + \ (d, X f - \m 2 X 2 + g&X ■ (1-1) 

In order to compare with standard literature (e.g. Chap. 10 in Itzykson & Zuber 
|35| ) it is advisable to switch to charged particles and to allow for different 

frequently only the exact solution of the ladder BSe for massless exchange particles is called 
the "Wick-Cutkosky model". Here we use this designation in a more general sense for a class 
of scalar models described by Lagrangians of the form 11.111.21 . 
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masses of the constituents. Then one needs two different types of nucleon fields 
and therefore the model which we are going to use is described by the following 
Lagrangian 

i.„ i 



dM 2 - M? \$ 



+ 2^) 2 - 2 m V + ^ 9i\*i\ 2 X- (1-2) 

l _ i=i 

The coupling constants g\ are now twice as large as in the neutral case as can 
be seen by collecting all pieces quadratic in the nucleon fields. Previously we 
have considered |36| a bound nucleon-antinucleon system in the neutral model 
described by the Lagrangian while now we will investigate a bound two- 
particle state described by the Lagrangian (|1.2|) . Both models are, of course, a 
far cry from reality and in the present treatment basically equivalent. Of course, 
there are differences from the fact that in nucleon loops two types of nucleons 
run around, i.e. that the determinant factor appears twice. However, since we 
will work in the "quenched" approximation in the following this doesn't matter. 
Note that also the meson properties are unchanged in this approximation. For 
simplicity the present work will be restricted mostly to the equal-mass case 

Mi = M 2 = M (1.3) 

where Mq denotes the bare mass of the nucleon. First results from this new 
approach have been presented in ref. Here we will give a full account of 
this formalism, test it in the nonrelativistic case, present improved numerical 
calculations and explore special cases, in particular when self-energy and vertex 
effects are switched off but crossed-ladder diagrams are retained. 

Finally it is worthwhile to recall the historical and conceptual roots of the 
worldline variational approach: for a single scalar particle it may be considered 
as a relativistic version of the polaron problem, and hence the present work is 
close in spirit (but not in the details) to investigating a bipolaron, the possible 
bound state of two electrons in an ionic crystal. This occurs when the attractive 
electron-phonon coupling overcomes the repulsive Coulomb interaction and has 
been the topic of many theoretical papers j33 IHEj as a proposed mechanism for 
high-temperature superconductivity. 

This paper is organized as follows: in the next Sect, we briefly overview 
the well-known nonrelativistic results for the present model. We then study a 
particular current-current correlator (or polarization propagator) whose poles 
will give us the relativistic bound-state energies. Sect. 4 describes the variational 
calculation of that object in the worldline formalism, including the retarded trial 
action which we use. In the next two Sects, the variational equations for the 
retardation functions which enter our trial action are derived and numerically 
solved. The last Sect, contains our conclusions and the outlook for further work. 
Most of the technical details such as the calculation of the various averages with 
the trial action or numerical details are collected in several appendices. 

2 Nonrelativistic Results 



A useful guide for the results one may expect for the two-body relativistic bound- 
state problem is provided by the nonrelativistic limit. Although this is quite 
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standard we briefly recall the essential facts and point to the differences in a 
relativistic treatment. It is well-known (and re-derived by path integral methods 
in appendix IA.1J) that in the nonrelativistic limit one obtains a Schrodinger 
equation for the two particles interacting via an attractive Yukawa potential 

y( Xl -X 2 ) = - 1 ^-'Hxi-x 2 | ( 2-1 ) 

|X1 — x 2 | 



where 



= 9192 (2 2) 

16vrMiM 2 1 ' ' 



is the usual (dimensionless) coupling constant 2 . One of the simplifications of the 
nonrelativistic treatment is that center-of-mass and relative motions decouple. 
For massless mesons (m = 0) the Yukawa potential in Eq. (|2.1|) becomes an 
attractive Coulombic potential with internal energy levels 



nonrelat. 



a 2 



2(„ + l)* M »" — M.2... (2.3) 



where 

M red = MM (2.4) 
Mi + M 2 

is the reduced mass of the system. If m ^ then first-order perturbation theory 
gives 

e = -i a 2 M Ted [ 1 - 45 + 6<5 2 - 8<5 3 + O (5 4 ) ] (2.5) 

where 

* = -^r- (2-6) 

2aM red v ; 

is proportional to the ratio of Bohr radius to meson Compton wavelength. This 
shows that the finite range of the Yukawa potential, as expected, decreases the 
binding energy. Indeed, numerical calculations [HHl demonstrate that an attrac- 
tive Yukawa potential only develops bound states if 

r 1 = > 1.67981. (2.7) 

m 

Thus for large enough coupling constants and/or small meson mass the WC 
model will most probably - inasmuch as the nonrelativistic approximation is an 
useful guideline - have bound states. Note that nonrelativistically the masses are 
always constants which are not changed by the interaction whereas relativistically 
the binding energy subtracts from the "weight" of the system. If the system has 
total momentum q then its total energy is 

Q 2 

M l +M 2 + 2(Mi + Ma) V(M X +M 2 + e ) 2 + q 2 (2.8) 



2 Sometimes, e.g. in ref. I35| , denoted by A. Taking into account the different coupling terms 
for the case of neutral particles the present definition is the same as in refs. |24| - |29| . 
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and with eo < the rest mass of a bound state is less than the sum of its 
constituent masses. Of course, the separation of center-of-mass and relative mo- 
tion is no longer valid in the relativistic case, nor is the concept of a reduced 
mass. However, our approach of calculating the poles of suitable Green functions 
does not require these, and indeed automatically gives the correct nonrelativistic 
center-of-mass energy q 2 /(2(Mi + M2)) as shown in appendix IA.2I 



3 The Current-Current Correlator 



"In the relativistic approach, bound states and resonances are identified by the 
occurence of poles in Green functions." (p. 481 in ref. |35j). Therefore in general 
one has to investigate a 4-point function of the form 



G^(x 1 ,x 2 ,x 3 ,x 4: ) 



(3.1) 



with no external mesons and i,j = 1,2. However, for the purpose of getting 
bound-state energies this object contains too much information and therefore it 
is advisable to consider the Fourier transform of a special 4-point function where 
arguments are pairwise identical. In the context of the model described by the 
Lagrangian (|1.2|) this could be 



%(g) : 



d x e 



iqx 



T 



^j(x)h(x) &I(0)$j.(0) 



°) 

/ connected 



(3.2) 



Here we already have used translational invariance to set the second argument to 
zero. Equation (|3.2j) is a version of the polarization propagator in the language 
of many-body physics |40| or a 2-point correlator for the "current" operator 



it/ 



it, 



(3.3) 



It is also (up to a sign and different currents) essentially the quantity considered 
by Shifman et al. in their work on QCD sum rules |3J to obtain hadronic masses 
from quarks and gluons. Note that we have taken a particular scalar current and 
that we are free to choose different ones: for example, ^ i (x)d IJi ^j(x) may be used 
to study bound states with angular momentum 1. Due to covariance the present 
polarization propagator can only depend on invariants, i.e. H-ij(q) = Uij(q 2 ). We 
will concentrate on the case i = 1, j = 2 in which the operator (|3.3j) creates a 
nucleon of type 1 and one of type 2 out of the vacuum. 

To see what the structure of the polarization propagator is we now look at 
its spectral representation. This is done in the standard fashion by inserting a 
complete set of states with four momenta P n 



1 = \Pn)(Pn\ ; Pl = M 2 n 



(3.4) 



into Eq. (|3,2[) . Using translational invariance one can perform the x-integration 
and obtains 



n^(g 2 



(2-) 3 E 



1 



Ml + i0 



q -p n ){ qo +p° n ) 



p 



<%■«» 
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showing that the polarization propagator indeed develops poles at q 2 = M 2 . 
Here M n is the mass of the intermediate state which can be reached from the 
ground state (the vacuum) by application of either the current Cjj(0) or by 

CV/(0). Using the standard expansion of the field operators in terms of creation 
and annihilation operators for the free-field case, one easily sees that for it creates 
a two-particle state (consisting either of two nucleons or two antinucleons of type 
i and j) as intermediate state when acting on the vacuum. 



P, 



C«(0) 



(3.5) 



3.1 Reduction to a Path integral Over Meson Fields 



With no meson source the generating functional for our model is denoted by Z' 
and reads 



Z'[J*,J] := JllV^V^Vx expji J d 4 



(3.6) 



One may integrate out the nucleon fields to obtain 

2 / 



Z'[J*,J] = const. / T>x e 



JSobc] 



m 



i 



x exp 



\\det Oi( X ) 
i I d 4 xd 4 yJ*{x)(x\Or\x)\y)Hy) 



where 



So\x] = 2/ d A x [(d X ) 2 -m 2 x 2 ] 
is the free meson action and 

Oi{x) = -d 2 -M 2 + g lX 



(3.7) 



(3.8) 



(3.9) 



is the operator for the quadratic part of the nucleon action. Note that charged 
nucleon fields lead to a power of —1 of the determinant instead of the usual 
—1/2 for neutral fields. However, in the quenched approximation where this 
determinant is neglected this makes no difference. 
The functional differentiation for 



d A x e iq - x 



<5 4 \nZ' 



5J j (x)5J*(x)5J i (0) 6JU0) 



(3.10) 



j=j*=o 



is most easily done by employing a cumulant expansion in terms of the sources 
in Eq. (|3.7[) . This gives 

InZ' = O{J ,J 2 ) + ^f- z2jl(f d ^) J*k(^)Jk(^)Ji(x3)Jl(x'd 

xR kl {x 1 ,x 2 ;x 3 ,x 4 ) + 0(J 6 ) (3.11) 
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where 



R kl (x,y;x',y') = {{x \0^( X )\ v) \OT\x)\ l/» 

-{{x\0^{ X )\y)){{x>\Oi\ X )\y'))- (3.12) 
Here we have defined the following average over the meson fields 



/ 



T>x exp 



iS [ X ] - ^TrlnOib) 



x J V X A(x) exp 



iS \x] -^TrlnO^x) 



(3.13) 



The functional differentiations can now be performed easily with the result 
Uijiq) = i j d A xe iq ' x Rij(x,0;x,0) 

i [ d A xe iqx [«x|Or 1 ( x )|x = 0> (x \oj\x)\ x = o)> 



-{{x\Or\x)\x = 0)){(x Oj\ X ) x = 0)> 



(3.14) 



This describes the propagation of a nucleon (or an antinucleon) of type i and one 
of type j from the space-time point to x. In between they emit and absorb all 
types of mesons which is represented by the functional integral over the meson 
field x- 

The last term Eq. (|3.14l) subtracts the unconnected pieces and consists just 
of the product of usual propagators. Since this term does not contain any poles 
besides the usual one-particle poles we can drop it when we search for additional 
bound-state poles in the polarization propagator. Taking i, j = 1, 2 we therefore 
will investigate the simpler form 



n(g) := = i j d 4 xe i ^((x\O^( X )\x = 0) (x \0^(x)\x = 0)> 
3.2 Worldline Description of the Polarization Propagator 



(3.15) 



We now derive the worldline formulation for H(q) given in Eq. (|3.15j) still for 
unequal masses. As in the case of the 2-point function this is only possible in the 
quenched approximation, i.e. by neglecting the determinant in the average (|3.13|) 



quenched 



J £>X exp(iS [x] 



£>X^(x)exp( iS [x\ 



(3.16) 



Only then can one perform the path integral over the meson field X after using 
the Schwinger representation 
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s2 



y 



exp 



a; 



x(T)=y 



x(0)=x 



T>x exp < i / dt 



X H Y X 

2 2k 



toice, i.e. for 6oi/i propagators. In this way one obtains 



U(q) = i J d 4 xe iq x £ 



dT\ dT2 
N2 ex P 



n 



2k 

2 ,T, 



M{T\ + M|T 2 



o 



exp < i 

- \J Xl (o)=o J { i= i 



(if; 



«0 -2 

— ^ 

2 * 



i=l 



Writing the argument of the last exponential as % J d 4 zb(z)x{ z ) with 



b(z) 



7T- Y] 9i / dtS(z- Xi(t)) 
2k o ~( Jo 



(3.18) 



(3.19) 



(3.20) 



the required functional average over the meson field can now be performed since 
it is just a gaussian integral 



Vxe 



iSo[x]+i(b,x) 



T>x exp 
const, exp 



(x, (~d 2 - m 2 ) X ) +i(b,x) 



-d 2 



m 



(3.21) 



Hence 



with 



(exp [i(6,x)]> = : exp( iS int [xi,x 2 



(3.22) 



-Sint [XI,X 2 ] 



i% £ 9i9j 

2 ,T, 



dt 



dt' (xi(t) 



-d 2 



V / 'dt [ 3 dt> [ ^HH'i(*)-';% ) 
•-fa 8^0 7o 7 (2vr) 4 p 2 _ m 2 + i0 ^ ; 



3 ko > reparametrizes the proper time and can be considered as "mass" of the equivalent 
quantum-mechanical particle. In numerical applications we will set this parameter (more pre- 
cisely: its euclidean counterpart) to 1 in which case all proper times have mass dimension —2 
as can be seen from Eq. (13.1711 . Another convenient choice is ko = M so that the proper time 
has the same dimension as the ordinary time to which it reduces in the nonrelativistic limit 

ED 
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The polarization propagator is therefore given by the double worldline path 
integral ^U] 

dT\ (IT2 



IL(q) 



where 



d 4 xe iq - x 







(2iK ) 



cxp 



-— (Mfc + M%T 2 ) 



n 



l \Jxi{0)=0 



VxA i So[xj] + iSintjx!, x 2 ] >(3.24) 



S [xi 



1,2 



(3.25) 



is the standard free action for each particle. 





Figure 1. Perturbative expansion of the correlator in the worldline formulation of Eqs. 13.231 
- 13.251 . Full (dashed) lines represent nucleons (mesons), filled dots the meson- nucleon vertices 
and circled dots the spacetime points x — and x where the nucleon lines start and end, 
respectively. 



The perturbative expansion of the polarization propagator is shown in Fig. 
n Note that the worldlines of the particles are parametrized by their proper 
times over which one has to integrate; thus crossed and ladder-type diagrams 
are included on an equal footing as are self-energy and vertex corrections. 

If Eq. (|3.23|) is split into terms with i = j and i 7^ j one sees that the 
former generate the self-energies of each particle, while the latter describe the 
interaction between the nucleons by exchange of (any number of) mesons. The 
vertex corrections come automatically due to different values of the proper times; 
for example, if one self-energy meson is already "in the air" when another one is 
emitted to the other particle. A typical diagram which is contained in the com- 
pact expression (|3.24j) is depicted in Fig. [2^, . A diagram which is omitted due 
to the quenched approximation (neglect of pair production) is shown in Fig. ^p. 
We normalize the path integrals by dividing and multiplying with 




(3.26) 
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which is easily obtained from 

-x(T)=x 

Vxe iSo[x] 

x(0)=0 



cxp { iT^— ) j - (J 
1 2ac 



d A p ( p 2 T 
exp — ip ■ x + % 



Thus 

n(g) 

where 



dTx dT 2 
16vr 2 7 (Ti + T 2 )< 



exp 



(2vr) 4 

g 2 TiT; 



2k 



(3.27) 



2k Ti + T 2 2kq 
J V(xi,x 2 ) exp(z£) 
/ V(x 1 ,x 2 ) exp(i,So) 



(MfTi + M|T 2 ) 



(3.28) 



Z>(xi,x 2 ) 



il(Ti)=a: 



xi(0)=0 



Z2(X 2 )=a; 



aj 2 (0)=0 



5 = 



x + S , So 



• a; + S [xi 
i=l 



Vx 2 (3.29) 



(3.30) 



With no interaction the ratio of path integrals in Eq. I)3.28|) is unity and one 
obtains the free polarization propagator. Without performing the proper-time 
integrals in Eq. (|3.28|) one sees in this case that the exponential vanishes, i.e. 
there is no damping anymore if the ratio r = T 2 jT\ of proper times fulfills the 
quadratic equation 

r 2 M| + r (Mf + M\ - q 2 ) M[ = . (3.31) 
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This only has positive, real solutions if 

q 2 > q 2 h = (Mi + M 2 f . (3.32) 

Hence there are (branch-point) singularities in the free polarization propagator 
starting from values of the external momentum q as given in Eq. 1)3.32 J) : this 
is just the familiar cut in the polarization propagator which describes the two 
particles in the continuum. It is a cut because we integrate over two proper times, 
or equivalently over all the poles generated by one proper-time integration (say 
T\) depending on r. In the interacting polarization propagator this cut also exists 
although it is modified in strength but not in position by the interactions between 
the unbound nucleon 1 and 2. However, for the bound state problem we are 
looking for an additional pole below qL. As in the case of the two-point function 
we expect that such a pole may be generated by the undamped integration of 
an exponential over one proper time (combination). 



4 Variational Calculation for the Worldline Polarization Propagator 



We now apply the Feynman-Jensen variational principle to the last factor in Eq. 
(|3.28j) . i.e. approximate 

J V exp(zS') J V exp(iSt) J V exp(iSt) exp(i(S — St)) 



J V exp(iSo) / V exp(i5o) 
^ J V exp(iSt) 
J V exp(iS ) 

In the last line the average is defined as 



• exp 



/ V exp(iS t ) 



i(S - St 



(4.1) 



S-St 



JV(xi,x 2 ) [S-St] exp(iS'i) 



/ f>(xi,x 2 ) exp(iS t 



(4.2) 



It is essential for our approach that the RHS of Eq. 1)4.1)) not only yields an 
approximation but is stationary w.r.t. arbitrary variations of the trial action St- 
Because the path integrals have to be solvable we are restricted to quadratic, 
but retarded (i.e. non-local) trial actions. 



4.1 The Trial Action 

In the following we will use the Fourier path integral method as employed orig- 
inally in refs. |24l 12*5] 4 . Although there exist more general and elegant treat- 
ments |29) . the expansion in Fourier modes has the advantage of being simple 
and straightforward. 

It is convenient to rescale the times in the free and the interacting part of 
the action 

t = T iT , t' = Tjt' (4.3) 
4 In the following these references will be denoted by (I) and (II), respectively. 
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with r, t' G [0, 1]. The paths fulfilling the required boundary conditions can then 
be written as 



Xi T = XT 



fc=i 



Li „(*) 



A:7r 



(4.4) 



where, of course, the Fourier coefficients a y k ' are 4-dimensional vectors. The free 
action becomes 

So - 



K0 \ - 

2 ^ 

k=l 



,0 2 



i=l 



5 = q ■ x + S . 



This serves as a guide for a quadratic trial action which we choose as 



(4.5) 



St 



Xq 



x 



Ko 

2 



E 



l k 



,0 2 



fc=l 



(4.6) 



fe=i 



where the last term accounts for the direct coupling of the two worldlines. Being 
diagonal in Fourier space this is not the most general quadratic action but in 
all known cases (polaron, WC model, QED) possible non-diagonal pieces do not 
contribute in the limit of infinite proper time(s). However, as shown in ref. |29| . 
additional Lorentz structures built from external momenta (here the momentum 
transfer q) can lead to considerable improvements. Here we will not include such 
refinements and stick to the simple trial action (|4.6|) , As shown in appendix IA.2I 
for the nonrelativistic case, such a trial action leads to the correct center-of- 
mass energy and to an internal ground-state energy of the system exactly as if 
evaluated with a variational trial wave function of gaussian type. It may also 
be noted that in the relativistic case the ansatz (J4.6D generalizes the one for a 
single nucleon, which showed the instability of the WC model in the one-particle 
sector. 



4.2 The Pole of the Variational Polarization Propagator and Mano's Equation 



Since the trial action is at most quadratic in the dynamical variables, all path 
integrals needed for Feynman's variational principle can be performed. The cal- 
culation of these averages is outlined in appendix B and gives the result iJTTmJl 



n var (g) 



i 



dTi dT 2 



+ 



167T 2 

TtT 2 
Ti +T 2 



exp< 



i 

2k 



-(MfTi + M|T 2 ) 



(Ti+T 2 )2 

q 2 (2A - A 2 ) - (Ti + T 2 ) {Slu + V) (Ti, T 2 ) 



(4.7) 



where f2i2 and V are defined in Eq. (|B,17|) and (jB.18|) . respectively. The original 
variational parameter A has been replaced by 



(4.8) 
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As in the case of the 2-point function (or the nonrelativistic polarization 
propagator) the variational expression simplifies considerably on the pole. How- 
ever, now we have two times T\,T2 and it is not clear a priori which one, or 
which combination should go to infinity in order to produce the (single) pole. 
For the equal-mass case (J1.3JI the solution simply is to introduce 



Ti+T 2 



Ti -T, 



2 ■ 



(4.9) 



Since TiT 2 /(Ti + T 2 ) = T/2 - s 2 /(8T) we have from Eq. PT1) 



n var (g) 



where 



64vr 2 



dT exp< 



4Kn 



-4M 2 +q 2 (2A-A 2 ) + ^ln 5 (g 2 ,T) 



(4.10) 



g(q 2 ,T) = ^ 



k2T 



-2T 



ds exp< 

I 2k 



q 2 (2\ - X 2 



-2 7' ((),, + I ) [T+ S -,T- S - 



(4.11) 



It is obvious that an additional pole develops if In g(q 2 ,T) has terms linear in T 
for large T. This happens if both Vt\2 and V are independent of T and s in the 
limit T — > oo. Then 



4ko 



T^oo 



lng(q 2 ,T) ™ -4 ( 12 + F ) - ^ + O 



4k hrT 2 
"1 T~ 



(4.12) 



because the remaining s-integral is convergent. 

The proof that Q%2, V become constants in the large T-limit follows along the 
same lines as in the 2-point case: for T — > oo we may assume that the coefficients 
Ak = A(kn/T) and that the sums over Fourier coefficients become integrals over 
E = kn/T 



k 



A'TT 

~T 



dkf\A 

In this way we obtain from Eq. l|B.17j) 



T 

TT 



dE f (A(E)) 



(4.13) 
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duo 
2ri 



dE 



]n\A.W(E)AW(E!)- B 2 (E) 
AW(E)+AW(E) 



+ 



A( 1 )(E)A( 2 )(E) - B 2 (E) 



(4.14) 



When the coupling function B{E) vanishes we find that ^1\2\b=o = ^i + ^2 ; i-e. 
the kinetic term reduces to the sum of usual kinetic terms for the self-energy of 
the individual nucleons (see ref. (I)). 
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To evaluate V from Eq. (|B.18|) we go back to the original (unsealed) times 
t\ = T\T,t2 = T 2 t' and note that for 



h+t 2 



a := ti — t2, £ : = => / dt\ I dt 2 . . . = da dH . . . 

2 Jo Jo J-T -Vl/2 

(4.15) 

The integral over the relative time a remains whereas the integral over the total 
time £ produces the needed factor T. Therefore we obtain 

v =t^ cy /(^r»-^+io exp {i &w - xp • ^ } ■ 

(4.16) 

This only works if the quantity p 2 j defined in Eq. (|B.15|) becomes independent 
of the total time £ in the limit T — > 00. Inspection of Eq. (jB.15j) together with 
the addition theorems of trigonometric functions shows that this is indeed the 
case: 

^(tuh-^Tz) = 4(a,S; S ,T) T -=3° ^(a) (4.17) 

with 

*W ^ ! f g ^ - B 2 (E) { (f ) 



+ (1 " <%) 



AW(i£) + A^(E) , . 
i— ^ - B(J5) cos (Ea) 



.(4.18) 



This is the generalization of the pseudotime /i 2 (a) for the self-energy of a single 
nucleon encountered in previous work. Even the diagonal (i = j)-part now de- 
pends on both profile functions as well as on the coupling function B{E). If the 
coupling vanishes then /i| reduces to the usual pseudotime for the i th particle. 
Similarly, the off-diagonal (i 7^ j)-part also depends on both profile functions 
AW ; ^4(2) . The mutual dependence reflects the change of properties of the one 
particle in the presence of the other ("medium effect") and the vertex correc- 
tions mentioned before. Note also the missing factor 2 in front of the last term in 
the exponential of Eq. Q4.16|) compared to the single particle (2-point function) 
case. This means that on the average each particle in the momentum loop of 
the polarization propagator has momentum q/2 and - since (q/2) 2 < M 2 in the 
bound state - is slightly off-shell. Finally, it is evident that /if- = p?^ . 

Collecting all terms in the exponent of Eq. (J4.1(Jj) which are linear in T we 
obtain (what we usually call) Mano's equation 

= -4M 2 + q 2 (2A - A 2 ) - 4 ( fi 12 + V ) . (4.19) 

determining a possible additional pole in the polarization propagator. In view of 
the previous remark it may also be written as 

|) 2 (2A - A 2 ) = M 2 + n 12 + V (4.20) 
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which very much resembles Mano's equation for the 2-point function. The missing 
factor 2 on the RHS in front of the kinetic and potential terms is contained in 
the definitions (|B.17|B.18|) . One may split up the interaction as 



V = : V U + 2V 12 



(4.21) 



with 



Via 



9i9j 
4k 



da 



dp 

(27) 1 



l 



m 2 + iO 



exp 



75— (//4'0) -\p-qcr) 

ZKq 



(4.22) 

Here we have used the fact that the pseudotime is even to restrict the cr- 
integration to the interval [0, 00]. One may call Va the self-energy part and 
V12 the binding part of the interaction, although these effects are coupled as 
mentioned before. 

Let us discuss Mano's Eq. (j4.2()j) qualitatively: without the kinetic and po- 
tential terms one would obtain A = 1 as a stationary solution (as in the non- 
relativistic case) and q 2 = q 2 h = AM 2 would determine the beginning of the 
continuum branch-cut. The effects of the interaction can best be estimated in 
the euclidean formulation (kq = Ike). Then for small coupling 

2 







12 — 



dttE 



In A® + 



1 



A® 



1 



+ 



1 



+ 



1 



AW AM \AW A( 2 ) 



1 



+ 0(B 4 



(4.23) 



which is positive since In A + 1/A — 1 > for A > and the second term 
adds for A® < 2. Thus the kinetic term tends to push q 2 above the threshold 
value. On the other hand, from Eq. (|4.16l) the potential term can be seen to be 
attractive: V < since dpo = idp^ and p 2 — m 2 = —{p 2 E + m 2 ). This is the usual 
competition between kinetic and potential energy and one may therefore expect 
a solution of Mano's equation with q 2 < 4Mp hys , i.e. a truly relativistic bound 
state for sufficient strong coupling. Since we know that the instability of the WC 
model prevents large coupling, the nonrelativistic criterion (|2.7|) does not apply 
anymore. 



4.3 Regularization and Renormalization 

As usual in quantum field theory divergences show up and have to be absorbed 
in physical, i.e. measured quantities (renormalization). Fortunately, the Wick- 
Cutkosky model is a super-renormalizable theory (Eq. (|2.2j) shows that the cou- 
pling constants Qi have the dimension of a mass) and in the quenched approx- 
imation only a mass renormalization is necessary. This has been done already 
for the 2-point function in refs. (I, II) and, indeed, the source of the divergence 
is again the self-energy of the particles. To be more specific we expect that only 
the diagonal pseudotime vanishes at small relative times: 

o (4.24) 
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In contrast, the non-diagonal part of the pseudotime stays constant (see ap- 
pendix C): 

Hl 2 (a) ^ const. (4.25) 

This means that for a — > 0,p — > oo the p-integrand for Vu (but not for V\2 
!) loses its exponential suppression which leads to a (large p or small a) UV- 
divergence. The solution is straightforward and simple (compared to a real, i.e. 
renormalizable theory with a dimensionless coupling constant): regularize the 
divergent expression, isolate the divergence and combine it with the (squared) 
bare mass of the particle. In ref. (I) we simply subtracted a particular term with 
the meson mass as a scale. 

Another possibility for regularization is a cut-off at small proper time (proper- 
time regularization) , i.e. 



da . 



da . 



l/A 2 



(4.26) 



It is understood that at the end of the calculation we let the cut-off A (with di- 
mension mass squared) tend to infinity. We then write the regularized expression 
for Vu as 



Vu 



9j 

4At 



oo 

l/A 2 
oo 



da 



1 



d,a 



l/A 2 



d 4 p 

(2vr) 4 p 2 - m 2 + iO 

d 4 p 1 
(2vr) 4 p 2 - m 2 + iO 



exp 



2k 



p 2 a 



exp 



2k 



P /%0) - Xp-qa) 



exp 



2kq 



p a 



V s "* + V ■ 

r ii • ' 11 



(4.27) 



The singular part can be evaluated easily: 



V 



smg 



A 

4Kn 



da 



7i/A 2 



9J 







32vT 2 



du' 
1 



d 4 p 
(2^) 4 

exp 



exp 



\2k 



[(p 2 — m 2 )u + p 2 cr] 



o 



2k 



-m 2 u' 



32vr 2 



e z Ei(z)(4.28) 



u' + l/A 2 

where z = im 2 /(2KoA 2 ) and -Ei(z) is the exponential integral (ref. [12], Chap. 
5). Using its expansion for small arguments (i.e. large cut-off A) one obtains 

prop, time reg 



ysirig 



32^2 



— ) +C 



(4.29) 



where 



C 



9? 



7£ + ln( 2^ 



(4.30) 



prop, time reg. 

is a constant specific for this regularization scheme. The divergent part plus the 
constant (or part of it) can be absorbed in the bare mass of each particle by 
defining the finite (intermediate) mass 



Ms 



16vr 2 



In I ^ ] +2C 



in- 



(4.31) 
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For gi = (72 = 2g, C = this is exactly the quantity 5 which has been calculated 
in refs. (I, II) by requiring a pole of the 2-point function at the physical mass M. 
The remaining, regularized interaction just provided a finite mass shift to the 
observed mass of the particle. Note that in F L reg we now can send the cut-off to 
infinity since the integrals are convergent by construction. In this way all traces 
of the cut-off have disappeared, or better are hidden in the masses Mj. As can 
be seen in Eq. 1)4.30)1 proper-time and other regularization methods with a cut- 
off should only be used in euclidean time (ko = ike) to make everything real. 
Other regularization schemes do not need that restriction: for example, after 
performing the cr-integration, one has in d = 4 — 2e dimensions 



V sm& 



9i ,.2e 



d d p 1 2ikq 



u Ze / £l : u (4 32) 

dim 4k J (2ir) d p 2 - m 2 + iO p 2 

where v is a mass parameter introduced for keeping the mass dimension of the 
coupling constant fixed in arbitrary dimensions. The p-integral is a standard one 
(see, e.g. appendix B in ref. [331) and one obtains 

^sing 9i { 4irv 2 \ e r(e) 



dim 32-7T 2 \ m 2 J 1 — € 

f ,/2 \ 1 

4.33) 



32vr 2 



- + 1 - 7E + ln(47r) + In ( ) + Ou ) 



which is independent of kq and real. As usual one may combine the bare mass 
with the 1/e-divergence leading to the minimal subtraction (MS) scheme or in- 
clude some of the constants (MS). In the Pauli-Villars regularization one stays 
in d = 4 dimensions but modifies the meson propagator 

I > I I . (4.34) 

p2 _ m 2 _j_ jQ p 2 — m 2 + iO p 2 — A 2 + iO 

This just amounts subtracting from Eq. 1)4.33 J) the same expression with m 2 — > 
A 2 . Hence one obtains the same result as in Eq. ()4.29j) with no additional con- 
stant. 

Whatever regularization one chooses, the divergent part of Va can be com- 
bined with the (squared) bare masses of the nucleons. Of course, there is the 
ambiguity of adding possible constants to the singular part of the interaction 
and subtracting them from the regular part, leading to different renormaliza- 
tion schemes. The scheme used previously did not have a constant term in Eq. 
()4.29|) and we will also choose this convention. However, it should be stressed 
that this only affects the relation between the unobservable bare mass and the 
intermediate mass and all physical results are independent of this choice. 



4.4 Special Cases 

Let us now consider the case where the two particles not only have the same 
mass but also the same coupling constant 

91 = 92 = 9' ■ (4.35) 

5 Note that in these references the intermediate mass has been called Mi. To avoid confusion 
with the bare mass of particle 1 we use Mi, etc. in the present work. 
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Their dynamics being the same we can assume 

A<V(E) = A^{E) = A(E) 

and Mano's Eq. reads 
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Here 

^12 

and 

Vij - 



M 2 = (|) (2A-A 2 )-0 12 -2(y i T s + y 12 ) 

2A(E) 



A 2 (E) - B 2 {E) 



(4.36) 



(4.37) 



(4.38) 



9_ 

4k, 



o Jo J (2vr 



) 4 p 2 — m 2 + iO 



exp 



2kq 
5\j exp 



(p 2 l4j(<r) - Xp-qa) 

(4.39) 



2k 



p (7 



for j = 1, 2 . The pseudotimes are related to the profile functions by 



4 f 00 I 

^12^) 



oo 

dE- 



A 2 {E) - B 2 (E) 
1 



sin 



(4.40) 



7T J E 2 A 2 (E) - B 2 (E) L 
It is very useful to introduce the combinations 

Afc(£) := A(E)±B(E) 
because then the kinetic term separates in two distinct pieces 

fiu = n[A_] + n[A+] . 

Here 



A(£) - B(E) cos(Ea) .(4.41) 



L J 2ivr J 



In A(E) + 



1 



1 



(4.42) 



(4.43) 



(4.44) 



is just the usual kinetic term encountered in the self-energy of a single nucleon 
(see (I)). The pseudotimes now become 



2 1 

2 t'°° 1 



1 



1 



sin 



A_(£7) A+(£7) 
sui 2 (£ct/2) cos 2 (£cj/2) 



Ea 
~2 



A-{E) 



+ 



A + {E) 



(4.45) 
(4.46) 



Note the appearance of the cosine function in the 'interaction pseudotime' fi 2 2 
which leads to a finite value at a = 0, in contrast to the 'self-energy pseudotime' 
which vanishes at that point. As explained before this requires a subtraction 
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in Vn but not in the interaction part V\i- We may perform the p-integration in 
the usual way by exponentiating the meson propagator (cf. Eq. ()4.28l0 . Then we 
obtain explicitly in euclidean time (kq = ike ) 



T/reg 



a 

2tt' 



da 



du 



1 



aXq/2 
m/in{a), —,u 



12 



Za 
' ~2tt 



M 



da 



Mi 2 (°") 



o 



-— e (m \fo, 0, u) 



due m/ii 2 (o-) 3 jr^,^ 



(4.47) 
(4.48) 



where 



qAq/2 
m nij{a), ; — - , u 



My 



exp 



2«e 



u , (Ag/2)V 2 

1 s-^-; — W 



eAu,a) 



(4.49) 



is the function used in refs. (I, II). Here g/2 is a shorthand notation for \fq* /2 
and the dimensionless coupling constant as defined in Eq. 1)2.2)1 is given by 



a 



9 



12 



(4.50) 



16irM 2 

In Eq. 1)4. 47|) we now explicitly use the renormalization scheme of ref. (I) , i.e. 
subtract a simple exponential in the proper-time integral which has the meson 
mass m as a scale. Subtraction at a different scale v (which, e.g. is needed for 
m = 0) would add the term 

■i 



— M 

2vr 

— M 

2vr 



00 i 

da — 
o Jo 

1 2 

m 



d;U 



exp 



a 



i/ 2 1 



u 



2ke u 

,2 , 



exp 



m 1 



2K£ U 



■ a 



m 



du In -tt = — M In 



2vr 



(4.51) 



on the RHS of Eq. 1)4.47)1 . This is exactly the same term as obtained in Eq. 
1)4.33)1 directly from dimensional regularization and amounts to chosing another 
constant C in Eq. 1)4.31)1 . Therefore one obtains a new intermediate mass 



M\v) = MI--M 2 lnf^ 

7T 



(4.52) 



which is independent of m but independent. Obviously, the scale v is arbitrary 
and the physical observables do not depend on it. Note also that O12 and V\j 
are independent of the reparametrization parameter ke- This is because profile 
functions and pseudotimes scale as 

A±(ke,E) = A±(k e E) =>- filj(KE,(r) = k e Hij (a / k e ) (4.53) 

where the last relation follows from Eqs. 1)4.45) I4*!4l)|) . Using the explicit forms of 
kinetic terms and potentials it is easily seen that these quantities do not depend 
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on the parameter ke- Mario's equation then implies that also the variational 
parameter A is reparametrization- invariant. 

In addition, in Eq. (|4.4H|) we have introduced an artificial strength factor Z 
which in the end has to be set to one but allows one to distinguish between the 
binding part and the radiative corrections. This is standard practice in QED 
bound-state calculations for atomic systems including hydrogen. In the present 
case, however, the factor Z has no physical meaning since it does not reflect the 
different coupling constant of one of the particles; otherwise we would have to 
include it in the corresponding self-energy of this particle. With this artificial 
parameter to switch on and off one now can study several special cases: 

(i) Z = : Since there is no dependence on fi 2 2 and /xfj is symmetric in 

A^, A + (see Eq. (|4.45|) ) both profile functions give the same contribution 
and Mano's Eq. (J4.37|) for the 2-body case becomes identical to the one for 
the self-energy of a single nucleon. 

(ii) a = 0, Za / : This corresponds to the case studied by practically all 
relativistic bound-state approaches (with the exception of atomic systems 
where one does perturbative QED calculations): neglect self-energy and 
vertex corrections. Naively one then expects no sign of an instability in the 
WC model but we will see that this is not the case. Since now M = M 
Mano's Eq. simplifies to 

M 2 = (|) 2 (2A - A 2 ) - Q[A-] - n[A+] -2V U . (4.54) 

Still, when varying this equation one includes all crossed exchange diagrams 
(in variational approximation with our quadratic trial action ...), not just 
iteration of the ladder diagrams. It is unclear how one can simulate the 
ladder Bethe-Salpeter approximation in the present worldline approach as 
it naturally contains all orderings of internal lines. Similar remarks apply 
to approximations where a 3-dimensional reduction has been performed 
jSJ) or one particle has been put on the mass-shell [7]. 

(iii) c — > oo : For the nonrelativistic limit one could re-introduce the velocity 
of light c and let it tend to infinity as in appendix IA.1I A simpler method 
is to perform the limit M — > oo in the unsubtracted expression (j4.22j) 
keeping the coupling constant a fixed. Setting q = (2M + e q ,q) one sees 
that the term p ■ q in the exponential forces po — > in all other terms. The 
Po-integration then produces a factor 

^ /A^A = ^o 2n6{ah (455) 
V 2k ) Ago 



i.e. an instantaneous pion exchange. Therefore 

_ 47rZJ- 1 aM 2 f d 3 p 1 
lj ~^ Xq J (2tt) 3 p 2 + m 2 6XP 



(4.56) 



Because /ifi(0) = the unregularized self-energy interaction V\\ now is 
just a (divergent) constant which changes Mq — ► M 2 but has no effect on 
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the dynamics of the particles (we anticipate A — ► 1 for M — ► oo). Writing 
/if 2 (0) = • V 1 ^ Mano's Eq. Q4.37JI then reads in the nonrelativistic limit 

M 2 = (M 2 + e q M-q 2 /4) 2 (2X - X 2 ) - - Q + - 2V 12 (lo) (4.57) 

from which we immediately deduce by variation that A-{E) = 1 as V12 
does not depend on it anymore. Hence J7_ = . Assuming V\2{oj) <C M 2 
one finds A — > 1 as expected and therefore for the total nonrelativistic 
energy 



AM + M 



Sl + + 2V 12 {u) . (4.58) 



Variation w.r.t. A + (E) shows that this profile function has the familiar 
(euclidean) form A+(E) = 1+lu 2 /E 2 , in agreement with the nonrelativistic 
variational calculation in appendix IA.2I (see Eq. ()A.50|) for M\ = M 2 )- 
Consequently, Mano's equation for the internal energy eo (with kq = ike) 



d f d 3 p 1 p 2 



60 = IM KEU; - 47TZa J (2^p^T^ eXP 



2ke w 



(4.59) 



is the correct nonrelativistic variational equation ()A.55|) for the Yukawa 
potential if we put d = 3. For the Coulomb case (m = 0) where all calcu- 
lations can be done analytically this gives the results 

k e uj = -f-(Za) 2 M 2 (4.60) 

d Z 7T 

However, there is no compelling reason to reduce the dimensionality to 
d = 3 by fiat and therefore we conclude that the nonrelativistic limit of 
the variational calculation is different from the variational calculation in 
the nonrelativistic limit which starts from the beginning with d = 3. This 
phenomenon can be traced back to the ansatz (|4,fi[) for the trial action 
which is covariant but too rigid to allow a different treatment of time and 
space- like dynamical variables needed for the correct nonrelativistic limit. 
An obvious remedy is to give the trial action more freedom, for exam- 
ple by allowing more general Lorentz structures for the profile functions 
A(E) -» A^(E) = A L (E) qf"q u /q 2 + A T (E) (g^ - q^q u /q 2 ) and similarly 
for B(E). Indeed, the kinetic term then splits into £li/d + {d — l)f2r/d 
(see Eq. (19) in ref. 29 ) which would provide the necessary dimensional 
reduction for the nonrelativistic limit without destroying the covariance 
of the whole approach. Actually, as shown in appendix C of ref. |29| . the 
variational principle demands such a general form if the profile functions 
are left completely free. 

a — > : Although our main interest is in the strong-coupling case a detailed 
investigation of the weak-coupling limit is worthwhile, for example to see 
whether the logarithmic terms of the ladder Bethe-Salpeter approximation 
(see, e.g. Eq. (10-78) in ref. j3S]) also appear in the present approach. 
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Since the nonrelativistic limit coincides with the weak-coupling limit we 
can expect that Eq. Q4.61|l with d = 4 is the leading term for m = . 
Indeed, as shown in appendix IF.ll the binding energy in the worldline 
variational approach has the expansion 



M/2 



(Za) 



7T 



7 a 

Z IT 



(Za) 4 



7T 



1 + 



(4.62) 



This shows an increased binding caused both by radiative effects (indicated 
by the a/V-terms) and by relativistic corrections (indicated by the (Za) 4 - 
term). Further discussion is postponed until Sect. 16.21 

5 Variational Equations 

Here we derive the variational equations for the A— parameter and profile func- 
tions A±{E) from Mano's equation ()4.37j) . As this is straightforward we just 
quote the results: assuming that Ke does not depend on A variation w.r.t. A 
gives 6 
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du u ej(u, a) . (5.2) 



The explicit A— factor in the derivative can be combined with the LHS of Eq. 
(|5.1|) to obtain the variational Eq. for A in the form (cf. Eq. (139) in ref. [I]). 
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This is still a nonlinear equation but is expected to better converge under itera- 
tion as part of the interaction is already included. In addition, Eq. (|5,3[) explicitly 
shows that < A < 1 since the integrand is positive. 
Variation w.r.t. the profile functions A±{E) gives 
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The functional derivative of Vu w.r.t. the pseudotimes is easily found to be 
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One can choose a A-dependent reparametrization parameter to simplify the potential energy 
term V at the price of making the kinetic term Q12 A-dependent. Due to the virial theorem 
(see below) nothing new is obtained by such a procedure. 
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which should be compared with Eq. (140) in ref. [I]. With a suitable integration 
by parts the result can also be written as 
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and is particularly simple for massless pions: 
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Note that we have to deal with one-dimensional but highly non-linear integral 
equations. However, even without solving them they can be used to deduce the 
behaviour of A±(E), /xf - (a) for small and large values of its arguments. This 
is derived in appendix C. Here we just summarize the results: at large E both 
profile functions approach unity in the same way: 



The same holds for the pseudotimes in the limit a — * oo 



(5.9) 
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For E 1 , a — > there is an essential difference between the self-energy and the 
interaction part: whereas A-(E) approaches a constant value ^4-(0) , A + {E) 
diverges for small E: 
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Therefore this term cannot be considered as "small" but has to be kept con- 
sistently. This is, of course, what we expect for a bound state which cannot be 
reached by perturbation theory from the free solution. Similarly, the pseudotimes 
behave very differently for small a: ^f 1 (<r) -xras usual, but 
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tends to a constant. This different behaviour in the UV-region was crucial for 
the mass renormalization in the WC model. 

Given the solution of the variational equation it is possible to express the 
kinetic term O12 = f2_ + ^+ not only from its definition in terms of A±{E) but 
also from a "virial theorem" in which case it is given in terms of the pseudotimes. 
Indeed, following the derivation given in appendix E of ref. |29| one can derive 
straightforwardly that 
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Note that this relation only holds after variation and for the sum of the two 
kinetic terms. It is also possible to derive expressions for the individual terms 
but they are more involved and will not be considered here. The factor of "2" 
in front of Eq. 1)5. 13)) is due to our definition of /x^ - in Eqs. ([4.4514.46)) . whose 
normalization is reduced by one half (compared to the usual definition of the 
pseudotime) but where one has to sum over two profile functions. Indeed, if 
either V\\ = or V\2 = one obtains the standard virial theorem in the polaron 
approach (see ref. (II)). 

Another useful relation is obtained by differentiating Mano's equation with 
respect to some parameter, say the artificial coupling strength Z, and using the 
variational equations to greatly simplify the result. For this purpose we first note 
that the interaction potentials depend on the combination 

x:=\l^X^f (5.14) 

and, of course, functionally on the pseudotimes, i.e. on the profile functions. 
Introducing this combination as a variational parameter instead of A, Mano's 
equation therefore reads 

M 2 = qx - x 2 - ( tt+ + n_ ) - 2 ( V™ s + V 12 ) • (5.15) 

Note that Q± is a functional of the profile functions which, after solving the 
variational equations, are complicated functions of Z as are the parameter x and 
the momentum q. Thus by differentiation the RHS of Eq. (|5.15|) w.r.t. Z we 
obtain 

-Ef { ppj K + «- + ^ + v a) ] } 

At first sight this looks rather complicated, but due to the variational equations 
the terms in the curly brackets vanish identically and after reintroducing the 
parameter A we obtain the simple result 

-dz\V =A^' (5 ' 17) 

This is our version of the "Feynman-Hellmann theorem" [2] for the worldline 
approach. Note that the RHS originates from the explicit Z-dependence of V12 
and that the dependence of the variational solutions on that strength only enters 
implicitly. We use the Feynman-Hellmann theorem in appendix IF. ll to derive the 
first terms in the weak-coupling expansion of the binding energy. 



26 



Variational Worldline Approximation for the Relativistic Bound State 



6 Results and Discussion 

6.1 Numerical Solutions and Binding Energies 

The system of coupled integral equations (|5.3|) - (|5.5jl . together with Eqs. IJ5.7|) . 
Q4.45|) and ()4.46|1 can be solved by iteration using similar methods as in ref. (II). 
The basic idea is to use a grid of gaussian integration points Ei,aj on which 
profile functions and pseudotimes are evaluated and used in the corresponding 
E, cr-integrals. This is similar to the mesh methods employed in ref. [IS] for 
quantum-mechanical eigenvalue problems. The necessary modifications for the 
worldline bound-state case are discussed in appendix D. To obtain stable and 
reliable results it was crucial to understand and match the analytical solutions 
at small and large values of the variables (see appendix C) to the numerical 
outcome. 

After considerable effort we have developed a program which solves the 
variational bound-state equations numerically with sufficient accuracy. Setting 
Vu = reproduced the one-body results - a necessary but unfortunately rather 
weak check of its correctness. More stringent is the "virial check" , i.e. comparing 
the total kinetic term ()4.43j) evaluated from the profile functions with the one in 
Eq. 1)5.13(1 obtained from the pseudotimes. To find a bound-state solution one has 
to search for values q 2 < AM 2 where the LHS minus the RHS of Mano's equation 
(|4.37|) changes sign. This was done by applying the regula falsi once the possible 
solution had been bracketed. The values of the intermediate renormalized mass 
M have been recalculated at the beginning of the program by setting Za = 
and ^fq 2 = 2M. 




Figure Oh shows the behaviour of the profile functions A±(E) as functions 
of the variable E after a bound-state solution had been found in this way at 



\/q% = 1.99218M for a = 0.5. As expected from our analytic results the two 

profile functions have a very different behaviour for small E but it is also clearly 
seen that A_{E) differs appreciably from the one-body profile function obtained 
at the same coupling constant. This illustrates how binding affects the nucleon 
self-energies or - in other language - the in-medium effects. Similar changes can 
also be observed in Fig. Ob for the pseudotimes. 
Table ^ gives the result for the binding energy 



obtained for the mass ratio 



and a range of coupling constants a. The lowest value comes from the nonrela- 
tivistic criterion (|A.61|) which for the chosen parameters says that a > 0.4067 in 
order to get a bound state in the Yukawa potential (alternatively no solution of 
Mano's equation with q 2 < AM 2 can be found for too small a). An upper value 






(6.1) 



7 Our previous results |36) were obtained with masses M — 0.939 GeV and m = 0.14 GeV, 
i.e. m/M — 0.149095. In order to facilitate comparison with results in the literature we have 
repeated the calculations for this standard mass ratio. 
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Figure 3. (a) The profile functions A-(E) and A+(E), (b) the pseudotimes ^12(0") and ^^(cr) 
for the bound-state solution at a — 0.5 with reparametrization parameter ke = 1. For compar- 
ison the one-body results (subscript o) are also shown. 

is determined by the critical value 

«2t = 0.817 (6.3) 

above which no real solutions of Mano's equation in the one-body case are 
found 8 ; this is the variational sign of the instability of the WC model when 

8 This value is slightly larger than the one obtained in ref. (II) because of the slightly higher 
mass ratio. Similar values have been found with truncated Dyson-Schwinger equations |46ll3|. 
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Table 1. Results for the binding energy e = \/gf — 2M from the solution of the variational 
bound-state equations for m/M = 0.15 and different coupling constants. The value of the 
intermediate mass M of the single nucleon (renormalized at v — m) is also given. Calculations 
were performed with 3 x 72 gaussian points and the iteration was stopped when the relative 
deviation A (defined in Eq. (|D . 7| ) was less than 10 -3 or when it increased twice consecutively. 
In this case A m i n is given. The column "vir. check " gives the relative deviation (fiijf — fii2)/fii2- 
The width of the bound state above the critical coupling a y cr [ t — 0.542 is estimated from Eq. 

the quantum corrections are taken into account. 

However, already for a > 0.542 one observes that the iteration only converges 
for a certain number of iterations and then the global measure of deviation de- 
fined in Eq. (|D.7|) starts to rise again . In the one-body case this was a signal for 
the instability of the WC model and here also we can interpret this phenomenon 
as the impossibility of finding real solutions of the variational equations for cou- 
pling constants larger than 

«St = °- 542 • ( 6 - 4 ) 
Consequently the virial check deteriorates rapidly. We also observe that this 
(infamous) instability now occurs at lower values of the coupling constant. This 
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is not unreasonable as, in general (e.g. in a laser), transitions may be induced by 
the presence of other particles |37]: the higher the number of particles in a system 
the faster it decays. In ref. (II) we showed that above the critical coupling only 
complex solutions of Mano's equation were possible where the imaginary parts 
determined the width of the particle. Similarly, in Sect. I6.3l we will give a rough, 

(2) 

analytical estimate of at cv [ t and determine the width of the bound state above 
this critical coupling. 

6.2 Comparison with Other Work 

Of course, not only the possible width of the bound state due to the instability 
of the WC model is of interest but also the magnitude of the binding. Here our 
approach should be compared with the nonrelativistic variational method (using 
a non-retarded quadratic trial action, i.e. gaussian wavefunctions) and with the 
exact nonrelativistic results from solving the s-wave Schrodinger equation for a 
Yukawa potential in the reduced sytem. The results are listed in Table |2j 

In addition, we also give binding energies from the ladder Bethe-Salpeter 
equation (BSe). To be more precise we display the coupling constant belonging 
to a given < 2M calculated from Efimov's variational approximation [IB] 
to the ladder BSe. Comparison with some exact numerical solutions available in 
the literature (e.g. ref. [IH]) shows that this approximation with 2 variational 
parameters is surprisingly good. For example, for a binding energy of 0.01M 
ref. [50] reports a necessary coupling constant a = 0.5716 whereas Efimov's 
approximation requires a = 0.5730 9 . The 2-dimensional integrals needed in 
Eqs. (56), (58) of Efimov's paper were evaluated by standard numerical Gauss- 
Legendre integration and the minimum was found easily with the help of the 
CERN routine MINUIT. 

As can be seen in Fig. ^ the binding energies of the full variational approach 
including the self-energy and vertex corrections are much larger in magnitude 
than the results from the other approaches provided the coupling constant a is 
appreciably larger than the threshold value needed to support a bound state. 
Of course, in the present variational approximation this threshold value is too 
large due to the insufficient approximation of a Yukawa potential by a harmonic 
oscillator one. Also the comparison of nonrelativistic variational results with the 
exact values from the solution of the Schrodinger equation shows that quanti- 
tatively the variational approximation is not a very good one (for reasons that 
are discussed in appendix IA.2I after eq. IA.57|) . Nevertheless the field theoretical 
effects are clearly visible and produce much more binding. That this must be due 
to self-energy /vertex corrections can be inferred by comparison with Nieuwen- 
huis & Tjon's Monte-Carlo calculation [^2] which gives the two-body binding 
energy beyond the ladder approximation but still without self-energy and vertex 
corrections 10 . 

9 Actually, upon repeating the calculations listed in Table 3 of Efimov's paper for m/M = 0.1, 
better values were found in a few cases because the author had missed the true minimum |51|. 

Unfortunately, the results there are only given in the form of Fig. 1 from which values have to 
be read off in an approximate way: sfq*/M ~ 1.990±0.020, 1.981±0.013, 1.961±0.014, 1.921± 
0.018, 1.851 ±0.020, 1.770 ±0.016 for a = 0.4,0.5,0.6,0.7,0.8,0.9 (the abscissa in Fig. 1 is 4a). 
The BSe curve is in good agreement with the values obtained from Efimov's approximation. 
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Table 2. The binding energies e/M from the worldline variational approach with m/M = 0.15 
compared with those from the 3-dimensional nonrelativistic (NR) variational method, the exact 
solution of the nonrelativistic Schrodinger equation and (Efimov's variational approximation 
to) the ladder Bethe-Salpeter equation (lBSe). The last column gives the Monte-Carlo results 
(estimated errors of the last digits in parenthesis) of ref. (52| in which ladder and crossed-ladder 
graphs have been included ( "generalized lBSe" ) but self-energy and vertex corrections were still 
neglected. 
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Figure 4. The binding energy e/M of the two-body bound state for the equal- mass case with 
m/M = 0.15 as a function of the dimensionless coupling constant defined in Eq. (|4.50pi . The 
results from the worldline variational approach including self-energy and vertex corrections are 
compared with those from Efimov's variational approximation to the ladder Bethe-Salpeter 
(IBS) equation 48 . Also shown are the Monte-Carlo results (with errors) from ref. |52) in 
the "generalized ladder approximation" to the Bethe-Salpeter equation (glBS). The thick bars 
for the worldline results represent the width of the bound state above the critical coupling 
estimated from Eq. 16.251 . 

Having introduced the artificial strength parameter Z for the direct interac- 
tion it is, of course, possible to check the self-energy effects on the binding energy 
by considering the limit 

a — ► , Za — > a' (6-5) 

in which the self- interaction V\\ vanishes. For the single nucleon we then have 
M = M and the variational equations (J5.41 15.5j) simplify slightly because 

A-(E) +A + {E) 

Also the asymptotic behaviour of the profile functions changes to A±{E) — * 
l + u>lJ{2E 2 ). 

Fig. shows a comparison of variational binding energies with and without 
self-energy effects as functions of the coupling constant a! . It is seen that the 
binding energies without self-energy effects are much smaller especially near the 
threshold where the interaction is strong enough for a bound state. In this case 



2 + % • (6-6) 




Figure 5. Variational binding energies with and without inclusion of self-energy effects for 
m/M = 0.15. For comparison the results from the ladder Bethe-Salpeter equation (denoted by 
"IBS") are also shown. 



the particles move slowly, and without self-energy and vertex corrections the 
minimal coupling is therefore the same as in the nonrelativistic case. However, 
we have to correct for the fact that our trial action is too restricted to reduce 
to the correct 3-dimensional kinetic term (as discussed above) and therefore we 
have to use Eq. (|A~6"T1) with d = 4. This gives 



0.54207 . (6.7) 



Our numerical results obtained by solving the variational equations and collected 
in Table El are consistent with this minimal value, which is an additional, nice 
check on the correctness and stability of our numerical procedures. Note that the 
minimal coupling Q6.7|) for a bound state without self-energy effects numerically 
coincides with the value (|6.4j) found for the instability of the bound system with 
self-energy and vertex corrections. This points to an important role of the direct 
binding as a new "doorway" for the decay of the system. 

As noted before, the variational calculation "without" still contains (approxi- 
mately) all crossed diagrams and therefore leads to more binding than the ladder 
Bethe-Salpeter equation for sufficiently large a'. Note that up to a' = 1 no sign 
of an instability is found. However, for slightly higher coupling constants the 
variational calculation does not converge anymore and we find again a critical 
coupling 

1.041 . (6.8) 
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Table 3. Same as in Table □ but for a = 0, Za = a', i.e. without self-energy and vertex 
corrections which implies M = M. The width of the unstable state above the critical coupling 
o^rit — 1-041 is estimated from Eq. I|6.29|l . 



At first sight this is surprising as the instability of the WC model has been 
attributed to the self-energy effects which are absent in this limit and one would 
suspect a numerical instability. However, Table shows that the virial test 
(checking how well the variational equations are fulfilled numerically) is excel- 
lent up to the critical value. In addition, the very same analysis which explains 
the induced instability in the full model (with self-energy effects) leads to the 
prediction a' ~ 0.814 for zero meson mass (see next Sect.). We are thus led to 
the conclusion that the WC model is "unstable in nearly any case" except for 
rather special approximations like the ladder approximation. The Monte-Carlo 
calculations by Nieuwenhuis & Tjon |S2] may easily have missed this instability 
given their large error bars, systematic errors and the fact that they only consider 
coupling constants not larger than a' = 0.9. 

The increased binding due to the self-energy /vertex corrections is also evident 
in the weak-binding case where (for m = ) we have found in appendix IF. II 



M 

with an enhancement factor 



(Zaf (Za)' 

A var -| 7. h ■ . . 



7T 



(6.9) 
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K vav = 1 + -- + ... (6.10) 

Z 7T 
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Figure 6. One-loop corrections: (a) self-energy, (b) vertex. 



for the leading nonrelativistic term. That this is not an artefact of the variational 
approximation but also occurs in the full theory can be seen as follows: the one- 
loop vertex correction and the wavefunction renormalization for a free nucleon 
depicted in fig. give rise to an effective coupling constant 

g> — > g' eS = g' (l + ^ + ...) (6.11) 

as obtained in Eq. (58) of ref. [27j. As derived in appendix IF.2I this is also the 
exact one-loop result since the first-order variational calculation reduces to that 
in the weak-coupling limit. Thus we have a — > a(l + a/-K + . . .) 2 and we therefore 
expect that the exact bound-state energy of the equal-mass WC model has the 
weak-binding expansion 

M {Zaf 

e n = -—( — TTf ^exact + • ■ ■ 6.12 

exact 2 2(n + \y 

with an enhancement factor 

K exact = (1 + _ + ...) =l + 4- + ... . (6.13) 

The result ()6.1U|) from the variational worldline approach is in good agreement 
with this except that the numerical factor is 7/2 and not 4. Since the variational 
calculation underestimates the leading term - the nonrelativistic binding in units 
of the reduced mass is (Za) 2 /ir instead of the correct Coulomb value (Za) 2 /2 
- such a small quantitative discrepancy is of no concern. Note that the proce- 
dure used to determine the exact enhancement factor H6.13j) is equivalent to the 
matching procedure used in the application of effective field theories to bound- 
state problems |53| . Indeed, Eq. (2.2) in ref. [21] is exactly the nonrelativistic 
effective field theory for the massless WC model and in that language Eq. (|6.11|) 
determines the Wilson coefficient c\ to one-loop order. 

However, the enhancement observed above is at variance with Ji's claim |55j 
that inclusion of self-energy terms leads to less binding. This result was obtained 
in the light-cone formalism whose boosted states contain many equal-time Fock 
components. Therefore a direct comparison with equal-time methods is difficult 
although the physical results should be the same ... 
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The origin of the disagreement is not clear to us: it could be due to a miss- 
ing nucleon mass renormalization 11 or a shortcoming of Ji's restricted Tamm- 
Dancoff approximation which only retains (light-cone) states with at most one 
meson. In view of the fact that both variational and effective field theory predict 
more binding, serious doubts remain whether the light-cone calculation is cor- 
rect and/or complete. This assessment is also supported by the results from the 
dressed ladder approximation which give more binding compared to the calcu- 
lations with bare propagators (see Fig. 8 in ref. [3]). The observed enhancement 
of the effective coupling in the WC model is similar to the enhancement of the 
electron-photon coupling in QED: however, in the latter case only the magnetic 
coupling (the famous anomalous magnetic moment of the electron) is affected 
and not the electric charge which is conserved in a gauge theory. 

The C((Za) 4 )-terms in Eq. fUJl may be compared with the corresponding 
result from the Todorov equation, one of the relativistic quasi-potential equations 
which describe spinless binaries bound by a Coulomb potential. However, what 
is found in the textbooks (e.g. ref. Eqs. (4.197) and (4.206)) corresponds to 
scalar QED, not to our coupling of mesons ("scalar photons") to scalar particles. 
The latter case was evaluated by Brezin, Itzykson and Zinn-Justin who antici- 
pated Todorov's equation |58| in a study of the relativistic eikonal approximation. 
This approach includes relativistic recoil and "an approximate summation of the 
crossed-ladder Feynman diagrams" but "does not include radiative corrections of 
the self-energy type" (HHl- Using the (unnumbered) equation following Eq. (20) 
in that paper or Eq. (3.8) in ref. [SB] one obtains in the equal-mass case 



For n = this relativistic Balmer formula gives a coefficient 5/32 = 0.156 for 
the (Za) -term whereas the variational result (|6.9|) leads to 1/it 2 = 0.101. Thus 
both approaches predict an enhanced binding of the lowest-lying state also from 
relativistic effects. Although the variational calculation gives a smaller coefficient 
- as it does for the leading nonrelativistic C((Za) 2 )-term - it has the clear 
advantage of also containing self-energy and vertex corrections. 

Comparison with the weak-binding expansions shows that in the variational 
approach all calculated expansion coefficients have the same sign but are smaller 
in magnitude than in the exact case. Therefore it is tempting to deduce that the 
true binding energy is always below the variational bound-state energy as in the 
standard Rayleigh-Ritz variational principle of quantum mechanics. However, we 
have not succeeded in proving this conjecture rigorously. The main obstacle is 
that renormalization requires the intermediate mass M of a single nucleon as 
input in Mano's equation (|4.37|) and this quantity is the result of a variational 
calculation itself. Actually, as shown in Eq. (85) in ref. [I], M 2 is bounded from 

11 For a recent discussion of renormalization in light-front dynamics see ref. |5(i|. 



Todorov 




(6.14) 
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below and one may expect that a similar inequality also exists for the 2- 
body case. But this does not lead to a bound for the binding energy because 
this quantity is determined relative to the single-nucleon result. Note that this 
argument does not depend on whether the renormalization is finite or infinite: 
the minimum property of the variational calculation is also lost in the finite 
bipolaron problem |38j . 



6.3 Induced Instability 

In order to understand the numerical results it is also useful to have some rough 
analytical insight into the solutions of the variational equations. In particular, 
we would like to understand the mechanism of the induced instability. Such an 
insight is provided by the ansatz 

A-(E) = 1 (6.15) 

with A as the only free variational parameter. This was used in ref. (II) to de- 
termine the critical coupling beyond which no real solution of the variational 
equations exist anymore. However, in the two-body case the binding certainly 
has to be taken into account and the simplest choice is a harmonic oscillator-like 
form (ke = 1 ) 

A + (E) = 1 + ^. (6.16) 

This is also the nonrelativistic solution and has the correct form for E — > (see 
Eq. H5.11JI ). Hence there are two variational parameters to be determined from 
varying Mano's Eq. To allow an analytical treatment we furthermore assume 
m = and weak binding 

«1 (6.17) 



(q/2) 2 M 2 

together with the property that the relevant proper times are 0(1/ M 2 ) , i.e. 
loo" <C 1 . As shown in appendix E the critical coupling can then be worked out 
analytically and is given by 

vr (1 + Vl + 3z) 3 o 

For Z = 0, i.e. for the one-body case, this reduces to the value a cr i t = 7r/4 = 
0.7854 obtained in ref. (II). As a function of Z the critical coupling has the 
behaviour 



«crit ► J 



1 - -Z 2 + —Z A + ... 



(6.19) 



which shows that the critical coupling for the bound-state system is reduced 
compared to the one-body case. 

12 In the 1-body case the value of M served as a measure of the goodness of different parametriza- 
tions for the profile function. 
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Figure 7. (a) The critical coupling a cr i t and (b) the correction factor / cor r as function of 
the artificial strength Z of the binding interaction in the simplified treatment of the induced 
instability. 



Indeed fig. [T^t demonstrates that the reduction is monotonous with increas- 
ing Z . For Z = 1 one obtains 



"crit 



~ 0.4627 (6.21) 

Z=l 



which is in reasonable agreement with the value a cr it = 0.542 obtained from the 
numerical solution of the variational equations. It is also interesting that the 
low-Z expansion (|6.19j) does not converge for the value Z = 1, a phenomenon 
which has also been observed in bound-state QED calculations 60 . 

It is also possible to determine approximately the width of the bound-state 
for a > a cr it following the treatment in Sect. IV. A of ref. (II). To be specific we 
write 13 

q 2 = (M*- i^ , M* ~ 2M (6.22) 

and introduce 

C := A| ~ \(M-ij\ =: Coe-*. (6.23) 

A nonvanishing phase x signals a complex (approximate) solution of the varia- 
tional equations for coupling constants above the critical one. Near the critical 
coupling a cr it given in Eq. (|6.18|) the phase x(°0 is determined by the implicit 
equation (see appendix E) 



a = a crit 



[l + a 2X 2 + 0( X 4 )] (6.24) 



13 In relativistic calculations it is more customary to denote the square of the mass of an unsta- 
ble particle as M* 2 — iM*V but this does not make any significant difference for the rough 
analytical estimate which we will derive. 
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where «2 is given in Eq. (|E.26|) . This can be used to evaluate the width near the 
critical coupling for which one obtains 



r a ^ 1 2M ( ^— ^ ) f COII (Z) (6.25) 
with a correction factor 

ws) = ^(^^w+^f^y' 2 , , s ^ . (6 . 26) 

V ; (l + 2z + Vl + 3z) 3 / 2 

The limiting cases are 

(6.27) 

and fig. [TJo shows that this function is rising rapidly from / CO rr(0) = 1 to its 
asymptotic value / CO rr(°o) = 1.0607. At Z = 1 its value is /corr(l) = 1.08845. 
However, this enhancement of the width pales (and may be not even significant in 
the light of the drastic approximations which have been employed) in comparison 
with the factor 2 which distinguishes Eq. ()6.25|) from the one-body result. This 
factor 2 simply comes from the fact that the bound state has mass M* ~ 2M. 
Thus, the induced instability by the presence of the second particle not only 
shows up in the lower critical coupling constant but also in the width which is 
roughly doubled if one exceeds this coupling by the same relative amount as in 
the single-particle case. 

What happens in the limit (|6,5|) . i.e. for the case where the self-energy and 
vertex corrections are neglected ? This can be easily analyzed by writing Z = 
a' I a which goes to infinity when a tends to zero and cJ is kept fixed. With the 
help of Eq. (|6.20|) we immediately obtain 



2 (a — a cr it \ 



3/2 



9 tt a cr \t , 27tt 

°- sVsJC >=> a VI28 = °' 8141 < 6 - 28 > 

which is close to the numerically observed value. Similarly above the critical 
coupling the estimate (|6.25|) for the width becomes 

_ / a '_ a '.\3/ 2 

~ V2M , crit . (6.29) 

a=0 V "crit / 

Therefore the induced instability also shows up if self-energy and vertex correc- 
tions are switched off. 



7 Summary and Outlook 

We have investigated the relativistic binding problem of two particles with equal 
mass in the scalar Wick-Cutkosky model employing worldline variational meth- 
ods. Compared to the standard field-theoretical description the worldline ("par- 
ticle" ) representation of the system entails a huge reduction in the degrees of free- 
dom, leading to a quantum-mechanical path integral over the trajectories of the 
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heavy particles ( "nucleons" ) like in the famous polaron problem. While previous 
work in this approach concentrated on single-particle properties and processes 
this is the first application of the method to the inherently non-perturbative 
binding of two relativistic particles and can therefore be considered as a rela- 
tivistic analogue of the bipolaron problem: although the delicate balance between 
repulsive two-body interaction and attractive self-energy is missing, there is the 
additional problem of non-perturbative renormalization. Fortunately, the WC 
model is a super-renormalizable theory so that - in the quenched approximation 
- only an identical mass renormalization both in the 1-nucleon and the 2-nucleon 
sector is necessary. 

Whereas most bipolaron studies employed parametrized ansdtze for the vari- 
ational functions we used a quadratic trial action with free retardation functions 
which are determined by the variational principle itself. This variational princi- 
ple takes the form of a specific equation ("Mano's equation" for the two-body 
case) which is stationary under unrestricted variation of the variational functions 
and determines the lowest bound state in a fully covariant manner. No definition 
of center-of-mass coordinates, factorization of wave functions or choice of spe- 
cial frames is necessary. This was checked in the nonrelativistic case where the 
correct center-of-mass and internal energies were obtained from the variational 
calculation. In the relativistic case we have derived the (non-linear) variational 
equations, solved them numerically and studied the asymptotic behaviour of the 
variational functions for small and large values of their arguments. 

In this approach the binding energy is not determined by an eigenvalue equa- 
tion but by the specific value of the external variable q 2 < AM 2 where both sides 
of Mano's two-body equation are equal. We have obtained this value analytically 
for the special case of small coupling constant and zero mass of the exchanged 
particle ("pion") and numerically without any restriction below a critical cou- 
pling. In both cases we found more binding due to the effects of self-energy and 
vertex corrections which are included approximately in our calculation. This is in 
contrast to results from light-cone calculations [SSj which claim that self-energy 
effects are repulsive. However, a variety of results obtained by different tech- 
niques - dressed-ladder Bethe-Salpeter equation [3], effective field theories, our 
variational approach - all point to the opposite conclusion. We have identified 
the source of this increased binding as mainly due to the pion-nucleon coupling 
enhanced by vertex corrections, in a similar way as the anomalous magnetic mo- 
ment of the electron or the muon is enhanced by radiative corrections. Additional 
binding also comes from crossed-ladder contributions as shown in ref. [52]. This 
was confirmed by artificially switching off the radiative effects in our program 
and retaining only ladder and crossed-ladder contributions to the binding. 

However, whereas Quantum Electrodynamics has a stable ground state for 
not too strong fields the scalar Wick-Cutkosky model has not - a deficiency 
which too often is forgotten or neglected in bound-state calculations employing 
this model. Indeed, whereas previous worldline calculations in the one-particle 
sector already yielded a critical coupling above which no real solutions of the 
variational equations exist anymore, the two-particle bound state ceased to exist 
for even smaller coupling constants indicating an induced instability. Surprisingly 
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this instability also shows up when the self-energy and vertex corrections are ne- 
glected and only ladder and crossed-ladder diagrams are included. Although the 
variational equations in our approach are highly non-linear integral equations 
we were able to explain this phenomenon analytically by simple ansdtze which 
lead to solvable, algebraic equations. Actually, the ability to combine numerical 
results with simple analytical insights seems to be a feature which the worldline 
variational approach shares with ordinary Quantum Mechanics: the virial theo- 
rem which was used to control the accuracy of the numerical calculations and 
the Feynman-Hellmann theorem employed for the weak-coupling expansion of 
the binding energy are other examples. 

One of the main weaknesses of the present approach certainly is the use of 
a quadratic, oscillator-like trial action also for the direct interaction of the two 
particles which in the nonrelativistic limit reduces to a Yukawa potential. This 
limits the accuracy, in particular near the threshold where the coupling constant 
is strong enough to support a bound state. With the present (isotropic) trial ac- 
tion this deficiency is even worse since in the nonrelativistic (weak coupling) limit 
the kinetic term reduces to a 4-dimensional and not to a 3-dimensional theory. 
This could (and will) be easily improved by allowing anisotropic terms as was 
done in the one-particle sector in ref. thus giving the trial action more flexi- 
bility. The necessary modifications of Mano's equation are straightforward with 
a modest increase in numerical complexity due to the appearance of longitudinal 
and transverse components w.r.t. to the external momentum q. Second-order cor- 
rections to the variational result can also be calculated rather straightforwardly 
|61| and typically reduce deviations from known exact results by a factor of 3 - 
4. Still, one would like to obtain the correct quantum- mechanical binding energy 
for weak coupling also in this approach which could be possible following the 
work of Luttinger and Lu in the polaron case [H2]- At the expense of applying 
Jensen's inequality twice these authors have introduced a (variational) potential 
into the trial action and were able to reproduce the polaron energy for strong 
coupling where Feynman's quadratic trial action leads to the largest deviations. 

Further extensions may include the unequal-mass case, binding of more than 
two particles similar as in ref. |63| and inclusion of vacuum polarization effects, 
i.e. going beyond the quenched approximation (see ref. [HI] for attempts in this 
direction). Calculating bound states in QED would be another step towards 
approaching a realistic theory although the present-day perturbative techniques 
are sufficient for the small electromagnetic coupling constant. Spin degrees of 
freedom and gauge invariance are no obstacles for a worldline description [3*T] . 
As we were able to extract the anomalous mass dimension of the electron, i.e. to 
determine the singularity structure of renormalization constants in the worldline 
variational approach [HI], the necessary non-perturbative renormalization for a 
renormalizable theory like QED requires more work but should be feasible. It 
is clear, however, that the ultimate goal and challenge for a non-perturbative 
approach lies in the dynamics of quarks and gluons at low energies. Whereas color 
degrees of freedom can be included in a worldline formalism 65 a it is unclear, 
at present, how to treat the non-abelian - and therefore nonquadratic - gluon 
action in a gauge-invariant way. Auxiliary-field methods in which three- and 
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four-gluon self-interactions are eliminated 66 may be a possible way to proceed 
as can be shown for the simpler case of a $ 4 -theory jHTj. We think that the 
results obtained so far in the worldline variational approach, in particular the 
novel way of treating the relativistic bound-state problem presented in this work, 
are encouraging enough to make an effort in this direction. 

Acknowledgement. All diagrams were drawn by using JaxoDraw |ri8) . 



Appendix A: Nonrelativistic Binding 

A.l The Nonrelativistic Limit 

It has long been known that the exchange of a scalar particle (a meson with mass m) is 
approximately equivalent to an attractive Yukawa potential between the heavy nucleons with 
mass M. The usual derivation starts from the one-meson-exchange scattering amplitude in the 
static limit and equates it with the Born approximation of potential scattering (see e.g. ref. 
|69|. appendix 10). A more consistent way is to integrate out the meson field in the generating 
functional 



Z[J* ,J;j] = / J] 2>$ 4 V$: V X exp 



c 2 +^2(j*$i+$*Ji. 



-jx 



(A.l) 



for vanishing sources, keeping the velocity of light c (and for completeness also Planck's constant 
K) and then letting c — > oo in order to obtain the nonrelativistic limit (70) . The velocity of light 
(re) appears in the following places 



xo 
Mi 



ct 



9- 



c 2 dt 2 



(d'Alembert operator) (A. 2) 

(inverse Compton wavelength of nucleons) 
and it should be remembered that the action is the time integral over the Lagrange function 



S 



dt / dxC 



(A.3) 



entering as exp(iS/h) in the path integral. The standard shifted gaussian integral then leads to 
an effective action for the nucleons 



S, 
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which is a non-local $ 4 -theory. Making the ansatz |71| 
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and sending c — > oo 14 gives 

sur M — \ £ / d^; {x) (im + |U) - i E Iff 



-iMiC Z t/h 



(A.5) 



x y d 4 xd 4 y<j>*(x)(l>i(x) (a; | (-□ - m 2 c 2 /ft 2 ) x | y) <f)*(y)cj> (jj) . (A.6) 



14 Since the meson mass m may be small or zero we do not perform the limit in the meson 
propagator. 
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The first term is the kinetic term of a system of nonrelativistic particles (of two types) whereas 
the last one describes their interactions. In the limit c — > co the retardation of meson exchange 
(i.e. the dependence on po) is suppressed and one obtains an instantanous interaction 



-□ - m 2 c 2 /h 2 + iO 



J (27T 



exp[-jp- (x -y)/h] 



(2ttR) 4 pl/(c 2 h 2 ) - p 2 /h 2 - m 2 c 2 /h 2 + iO 
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.(A.7) 



ence, if the effective, nonrelativistic action is written in the standard many-body form for two- 
body potentials 



h 2 A 
2M~ 



0i(x,t) 



s cH [ft ,4>i] = J dt | j d 3 x <j>* (x, t) [ihdt + ■ 
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one finds that an attractive Yukawa potential 

h 4 gigj exp(-mc|x-y|/ft) 



Vij (> 



y) = 



l x -y| 



(A.9) 



with range h/(mc), the Compton wavelength of the exchanged meson, acts between the different 
types of nucleons. 

A. 2 Variational Approximation for the Nonrelativistic Polarization Propagator 

Since we make a variational approximation of the ground state energy of the relativistic system 
it is appropriate to assess the accuracy of our calculation in the nonrelativistic limit to which 
it should reduce for c — > co. The standard approach in Quantum Mechanics is, of course, 
the Rayleigh-Ritz variational principle which involves a trial wave function. Although it is 
obvious that a quadratic trial action corresponds to a gaussian trial wave function we will 
derive the corresponding approximation in close analogy to the field-theoretic (worldline) case. 
In particular, we do not use the simplification that the nonrelativistic Hamiltonian can be split 
up into a center-of-mass part and a relative part. Instead we will deal with the translationally 
invariant full system whose Hamiltonian reads 



* = ik + A +v ^-^ (A - 10) 

and investigate a quantity where one projects on a fixed (total) momentum q. This is just the 
nonrelativistic polarization propagator in d — 3 dimensions 



n(q,£) 
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= -i J d 3 xe-^ J dT (xi = x 2 = x |exp \iT [e - Hj] | xi = x 2 = OjJA.ll) 
which has poles at the total energy of the system 

n(q,£) = \M0)\ 2 {e 



_2(Mi + Mi) 

The time-evolution operator has the path integral representation 



(A.12) 



xi — X2 — x exp —iT [H — iO 



»)] 



Xl 



X 2 



= 
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V 3 xiD 3 X2 exp { iS [xi , x 2 ] } 
(A.13) 
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where the action is 



S , [xi,x 2 ] = 



dt 



Ml . 2 M 2 . 2 T/ , v 

— Xi + — x 2 - V (xi - x 2 ) 



(A.14) 



Note that both particles have one common time - in contrast to the relativistic case. This 
facilitates the limit T — -> oo even for the case of unequal masses Mi ^ M 2 . For vanishing 
interaction the free polarization propagator can be calculated easily 
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where M ro d is the reduced mass of the system defined in Eq. 12.41 . This is used for the normal- 
ization of the path integral in the general case. Thus we have 
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We now apply the Feynman- Jensen variational principle 14. 11 with a suitable (and manageable) 
trial action. For a proper choice we first look at the free action: with the Fourier expansion of 
the paths 

t ^ 2VT (i) . fk-nt\ 
a: ' sin 

I T J 



Xi(t) 



(A.20) 



the free action becomes 
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In the one-body case a rule of thumb is to "decorate" the free action with variational parameters 
to obtain a quadratic trial action which reduces to the free one if the interaction is switched 
off. In the two-body case we need an additional coupling term which should account for the 
interaction of the particles. Therefore we take as trial action 



-Aq-x + ^ 



,(1) „(2) 



(A.22) 



The Fourier coefficients Aq, as well as the parameter A are variational parameters. 

The relative minus sign between the A and B terms has been chosen for consistency with 
the relativistic ansatz l|4.6^ . Since the trial action is at most quadratic the averages can be 
calculated analytically: 
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Here 



V(j>) = J d 3 x V(x) e lp x 
is the Fourier transform of the potential between the two particles and the abbreviation 
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has been used. Thus we find 
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where the kinetic term is defined as 
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As usual the pole of the polarization propagator is determined by the limit T — > oo in which 
case the sums over Fourier coefficients become integrals over E = kn/T as in Eq. 14. 1311 . In 
addition, we may use the fact that only those terms in the exponent which develop a linear 
T-dependence contribute to the pole position. This means that the Ao-terms can be dropped 
and that we can replace 15 



sin (ET) = - [ 1 - cos( 2Et ) ] - 
in the argument of the function v(T,t) making it a constant: 
1 1 



(A.30) 



u(T,t) T -=5° v 



2tt MiM 2 J E 2 



;E 1 M 1 AW{E) + M 2 A™(E)-2VMZMZB(E) (a 



A( 1 )(E)A^(E)-B 2 (E) 



However, this constant is still a functional of the variational functions A^\E), B(E). Thus we 
obtain 



n 

v = 



3 
2iri 



dE 



In ( A w {E)A {: '> (E) - B*(E) 



(2), 



A m (E) + A (2) (E) 
AW(E)A&(E)-B 2 (E) 



d 3 p 
(2^)3 



V(p)e 



2(A.32) 



(A.33) 



Collecting the terms linear in T in Eq. 1A.28I and setting 

15 This can be seen, for example, by scaling t = T x and using well-known results for the asymp- 
totic limit of Fourier cosine transforms 1721. 
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Mano's equation reads 

Eo = 2[M f + Ma) (2A — A 2 ) + ( n + V ) =:ScM + eo . (A.35) 

Due to our metric (q 2 = q 2 — q 2 ) there is now a different sign between the center-of-mass (CM) 
and internal parts compared with the relativistic case. Variation with respect to A yields 

A = f (A.36) 

since O and V are independent of A. This is crucial both for obtaining the correct center-of-mass 
energy 

ScM = mrw) (A - 37) 

and for guaranteeing that the internal energy eo is independent of the total momentum of the 
system. 

To determine the optimal variational functions we now vary Mano's equation 1A.35I with 
respect to A W (E), A (2) (E) and B(E) 

jJw) {n + V) = °> sbW) {Q + V) = - (A - 38) 

The variation of the individual terms gives 

A (3-i) B 2 + A (3-i)2 



5Q 
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SO. 




Jb 


- «| 


8V 


C2 


SAW 


E 2 


8V 


C2 


SB 


E 2 



AWaW - B 2 (AW AW - B 2 f 



i = l,2 (A.39) 



+ ,\ „ I2 (A.40) 



AWAW-B 2 (AWAW-B 2 ) 

^mTI _ _ (yg>> + ^A^ + 2B) A^ 

(AW AW — B 2 ) 2 
-2 (a^A^ - B 2 ) + 2 [^A {1) + JW t A^ - 2B) B 
(AW AW -B 2 ) 2 



-i) 

-(A.41) 
(A.42) 



where we have omitted the argument E in A^'(E) and B(E); a = 3/ (2m) is a numerical 
constant whereas the constant 

1 1 f d 3 P 2f /t ,-i P 2 v 



C2 = ~~ , / jP^t p 2 V(p) e' lp " (A.43) 
depends on v, i.e. on the variational functions. We observe the relations 

2B SAW +A ( "IB = ^AWAW-Tb 2 - ( A - 44 ) 

2B JAW +A ^W = <»& V AWAW^B 2 S (A. 45) 
which (by adding suitable combinations of the variational equations) gives us 



i'VJi = /J.t', ( J^l + ^E 2 ] . (A.l.,1 

If this is inserted into one of the variational equations one finds after some algebra the simple 
solution 

B(E) = | ± (A.47) 



and hence 



""w-'+V^H- (A - 48) 
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Denning 



Mi + M 2 c 2 



VMiM 2 ci 

we finally have as solutions of the variational equations in real time 

■ 2 VMiMi uj 2 



A (i) {E) = 1 



Note that for equal masses 



1+ M i E 2 



M 3 



Mi + M 2 E 2 



A {1) {E) 



A (2) {E) =: = 1 



Mi=M 3 



so that A{E)-B{E) = 1 and +B(E) -■ 
for the quantity uj reads 

2 _ 1 Mi + M 2 
3 MiM 2 



1 u/_ 

2 £ 2 ' 

-uj 2 /E 2 . 



B{E) 
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2 £ 2 



(A.49) 



(A.50) 



(A.51) 



Mi=M 2 

With Eq. tQ31 the definition llA~49l 



/ 



d 3 p 
(2tt) 3 



pV(p)e- 



(A.52) 



This is an implicit equation for uj since v depends on ui. Indeed substituting the solutions <A,50t 
into Eq. 1A.3H one obtains 



2tt VMiM 2 



dE 



l y/Mi/M 2 + y%g 



E 2 



1 - w 2 /£ 2 + iO 
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4lu 



1 1 

M" + Ma 



(A.53) 



Here we had to specify how the singularity at E — ui has to be treated. This can be inferred from 
the original real-time path integral by noting that the profile functions need an infinitesimal 
positive imaginary part for convergence of the integral. This is equivalent to u — > uj' — iO, 
Feynman's prescription for the causal propagator. Note that the RHS of Eq. 1A.521 is positive 
for an attractive potential as anticipated. The implicit equation becomes more familiar if we 
consider Mano's equation before variation w.r.t. uj. The kinetic term 1A.32I is 



3 
2vn 



dE 



In 1 



E 2 



+ i0 + 



2 - uj 2 /E 2 



2 /E 2 + iO 



and therefore the internal energy in Mano's equation ljA.35^ is determined from 



3 



(27p 



V(p) exp 



4M rcd uj 



(A.54) 



(A.55) 



which is exactly the Rayleigh-Ritz variational principle with the normalized gaussian trial wave 
function 



/ -K/r \ 3 / 4 



exp 



M rod UJ 2 

_ — x 



(A.56) 



Indeed variation of Eq. IA. 5511 w.r.t. uj leads to Eq. IA. 5211 if Eq. 1A.53I is taken into account. 

Thus - as expected - in the nonrelativistic limit the Feynman- Jensen variational principle 
with a quadratic trial action is equivalent to the Rayleigh-Ritz variational principle with a 
gaussian trial wavefunction 16 . Apart from the mass-dependent factors the profile functions 
1A.5CH are recognized as the standard ones for a harmonic oscillator with frequency parameter 
uj . It is this variational parameter which has to be optimized for the specific potential V(x) in 
the quantum mechanical problem. 

For the Yukawa potential V(x) = — a exp(— mx)/x , x = |x| , a simple calculation gives 



to < e var {y) = MredQ 



d 1_ 

4 y 2 



_2_1 



+ 28 e" y erfc (Sy) 



(A.57) 



where y — ol^J M rc d Jul is the appropriate dimensionless variational parameter and S — 
m/(2aM rc d) was already defined in Eq. 12. 61 . erfc(a;) = 1 - erf(x) denotes the complementary 



16 As we have worked in real time where the Feynman- Jensen variational principle only ensures 
stationarity, we have missed the minimum property of the method. However, this is easily 
established in euclidean time where we can use the Feynman- Jensen inequality. 
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error function and we have explicitly introduced the dimension d(— 3) in the kinetic-energy 
term, d — 4 then also covers the weak-coupling limit of the relativistic variational calculation 
where - as discussed - the rigid, four-dimensional trial action leads to an enhanced kinetic term 
and therefore to less binding in the nonrelativistic limit. For small 5 one obtains 



M r ed a 



* +28 + . 
an 8 



(A.58) 



which means that the exact Coulomb result in Eq. 12.51 is missed by 15% for d — 3 and by 36% 
for d = 4. However, the first-order (for d = 4) and (for d = 3) also the second-order correction 
due to the finite meson mass come out correctly. At first sight this looks not very impressive for 
the variational calculation unless one realizes that a Yukawa potential is quite different from a 
trial harmonic oscillator potential. 
The variational energy becomes zero or positive when 

5 > 5 cvit = -±= (A.59) 

where so = Sy fulfills the transcendental equation 

s (1 + So) exp (so) erfc (so) - Sq = - . (A. 60) 

Numerically one finds so = 0.821324 and therefore <5 cr it = 1.10687/d. This means that the 
variational calculation predicts bound states for the Yukawa potential if 

2M rod f 2.71035 : d = 3 c1 , 
a ^T ~ { 3.61380 : d = 4 ^ AM ^ 

Due to the minimum property of the variational calculation this is always larger than the exact 
number $2.7) . 



Appendix B: Calculation of Averages 



To perform path integrals involving the trial action l|4.6^ we need the "master" integral in d 
dimensions 



hi 



d a x I V d a m V d a {2) e iSM 



V e 



(B.l) 



where 



2 oo 

c no , 2 K o v-> v-^ a(*) (*) 2 

Sm = p-x- —A x x - — 2^ 2^ A k a V + 

k=l i=l k=l 

This is practically the trial action with sources coupled to the modes which will be needed 
to evaluate the average of the interaction. With the (ubiquitous) gaussian integral 



i=l fe=l 



d d a exp \^~i^-Aa 2 + ib-a 
we get by completing the square 
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A d / 2 
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(B.4) 



where JT A / is a function of p, A X ,A^\ Bt and . 

We are now in a position to calculate the averages required in the Feynman- Jensen variational 
principle 14.11 . First 

f V exptiSt) 

= exp 



/ V exp(iSo) 



2k 



(JFt-jFo) 



(B.5) 
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To evaluate the numerator (i.e. Tt) one simply has to put p = Xq, A x = Aq(1/Ti + l/Ta) in Eq. 
l|B.4|l : for the denominator (i.e. To) one sets p = q, A x = 1/Ti + I/T2 and in both cases the 
sources are zero: = — 0. This gives 

/ V expjiSt) 
J V exp(iS ) 



exp 



1 

2k 



T1T2 



T 1 +T 2 V A 



X 2 



.(B.6) 



d = 4 is understood but kept general for pieces which contribute to the "kinetic" term in the 
variational functional. Next we obtain the average of So with respect to the trial action by 
differentiation 



So 



q a 2 00 

\ " \ 
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Similarly we may obtain 
1 d 
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Finally we calculate the average of the interaction part HM.2MJI 

(S^) t = ~E =i I dTdT 'J ^ jP-^ + iO < exp HP- (x a (r)-x,(r'))])- 



(B.9) 



Since 
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the average of the last exponential in Eq. HB.9t is given by 



exp 
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(B.ll) 



where 
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(B.12) 



With the help of the master integral l!B.4(l we obtain 



T\nt{a,f3) — T t 
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(B.13) 
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where 
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TiT 2 . 1,2 v-^ 2 
(r - r ) +2^ 



A)(T a + r 2 ) - ' ' • ^ (ftTr) 2 A m A m _ S 2 

■< <W r a fA< 2) 5 a i [sin(fcTrr) -sin(fc7rr' 



+ (1-<W) 



Putting everything together we now have 
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(B.18) 



Ti +T 2 

Equation 1B.16I still holds for all q, i.e. also away from the poles of the polarization propagator. 



Appendix C: Asymptotic Behaviour of Profile functions and Pseudo- 
times 

Here we derive the limiting behaviour of the solutions of the variational equations both for 
E, a — > and for E, a — > 00 . As the pseudotimes are Fourier cosine transforms of the (inverse) 
profile functions we expect, of course, that the large iJ-behaviour of the one function is related 
to the small er-behaviour of the other and vice versa. 

We first consider the limit E — > 00 in Eqs. (15.415.51 . Naively one would just replace the 
squared trigonometric functions by their average value 1/2 and obtain A± ~ > 1 + (consti + 
const 2 )/i? 2 where the constants are determined by dcrSVij/SfX 2 j(cr). However, for j = 1 
this integral does not exist as — * a for a — *■ and SVn/Sfiix diverges as 1/er 2 in that limit 
(see Eq. This is, of course the standard behaviour for the self-energy part and the way 

to obtain the proper asymptotic behaviour is well-known: introduce a factor a 2 /a 2 into the 
(j — integral and change to the variable x — Ea. We then obtain 

b^oc 2 a SVu f°°. sin 2 (x/2) , , aM 2 , „{\n 2 E\„, s 

k e E cr^o Sfif x (a) J x 2 4n E E \ E 2 / 
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where in the last step the (a — * 0)-limit of Eq. I|5,7^ has been used. It should be emphasized that 
this is exactly the same asymptotic behaviour as in the one-body (self-energy) case (see Eq. 
(143) in ref. [I]); the only new information here is that both profile functions become identical 
at asymptotic values of E. However, they have quite a different behaviour at small E: whereas 
A -(E) approaches a constant value 



A + (E) diverges for small E: 



1 + 



2k e 



da a' 2 



5Vl; 
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(C.3) 



as expected from the nonrelativistic limit. Consequently the pseudotimes also have quite differ- 
ent behaviour at small a. First, from Eq. 14.451 it is seen that /Ltfi has the standard self-energy 
behaviour 



— a f C 
T Jo 



dx 



\x/2) 



1 



A-(x/a) A+(x/a) 



dx ■ 



\x/2) 



For the low-cr behaviour of /j, 2 2 we may just set a = in the definition 14.461 to obtain 



dE 



E 2 A+(E) 



nh(o) 



(C.4) 
(C.5) 



where the integral exists due to the behaviour 1C.30 of A+(E) at small E. Thus for a — > the 
'interaction' pseudotime stays constant and nonzero as assumed. Note that this constant is not 
the inverse of the o; var defined in Eq. 1C.3II since A + (E) deviates from 1 + loI^/E 2 at E > 0. 
It is also possible to work out higher-order terms by considering 

1 1 



..2 , s 1 
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Because both profile functions have the same asymptotic behaviour ICH for large E the integral 
also exists for a = 0. Therefore we obtain after integration with the boundary condition 

A? 2 (o) = o 
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For the self-energy part the situation is different: 
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does not exist for a = as the integrand behaves like E for large E as can be seen from Eq. 
JUTj. This means that an expansion only in powers of a is not possible and that logarithmic 
terms appear in higher order. The leading term is obtained by realizing that the cosine provides 
an upper limit const. /a for the iJ-integration. Therefore 



p-ii(°) 





r- const, /fj 
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where ao is some undetermined constant. Integration with the appropriate boundary conditions 
gives 

(CIO) 



.0 , a AT 2 a 
-> a H a m — 

47T K E a\ 



where a\ is another constant to be determined at 0(a 2 ). In the framework of generalized func- 
tions one obtains the same result from the asymptotic behaviour IC.H of the profile functions in 
the definition 14.46H of [i\\(a) (see Table 1 from ref. |72| lr ). In turn this logarithmic behaviour 

17 Note the different definition of the Digamma function in this work compared to ref. |42| ! 
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of the pseudotime fi^ (a) induces the In 2 (E) / _E 2 -remainder in the asymptotic expansion l|C,10 
of the profile functions which is essential for the convergence of several integrals encountered 
above. 

Finally we consider the pseudotimes for large a by adding and subtracting the term 1/A_(0) 
from the inverse pseudotimes in the integrand. This gives 



MnO) = - 



dE- 



in 2 (Ea/2) 
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(C.ll 



where the last integral exists due to the proper behaviour of the integrand at E = and E — oo. 
Since the last line in Eq. IC.llll was obtained by replacing sin 2 (_Ea/2) asymptotically by 1/2 
it also holds for ^j.\ 2 where 1/A + (E) is weighted with cos 2 (_Bcr/2). 



Appendix D: Numerical Details 

Here we describe a few modifications and improvements of the numerical methods used in rcf. 
(II) in order to achieve stable and reliable results for the bound-state problem. As before we 
have used Gauss-Legendre integration rules with n g gaussian points and n e extra divisions to 
evaluate 

d v f(y) - \ w i f (va) . va = a + U - ^) h + \ x j > h = h —^- (D.i) 

i— 1 j=l \ / ~ 



where Xj,Wj are the usual abscissae and weights for Gauss-Legendre numerical integration in 
the interval [— 1, +1] (see Eq. 25.4.30 in ref. |42p. A major problem is the calculation of the 
pseudotimes from the profile functions at large E, a since the integrand sm 2 (Ea/2) leads to huge 
oscillations and loss of accuracy when the integrals are evaluated by mapping to a finite interval 
(e.g. by the transformation E = Eotanip,a = Eq tan</> , _Eo a suitable scale) and subsequent 
gaussian integration. This procedure was adopted in ref. (II) and found to be sufficient for the 
one-body case. Recently we have shown that one can eliminate the profile functions altogether 
and solve a nonlinear, delay-type equation for the pseudotime which bears a striking similarity 
to the classical Abraham-Lorentz equation |32| . This works perfectly well for the one-body case 
but would need a different strategy for the bound-state problem. Therefore for the present 
purposes we have decided to stay in the proven, conventional scheme but to take finite upper 
limits -E m ax, cr max for the integrations. The upper limit for E is chosen such that the asymptotic 
contribution 1C.H at E = -E ma x is less than 5 x A where A is the measure of deviation which 
should be reached in the iteration (see below). Thus for ke = 1 

(X M 2 

x ~ 4 5A ' ( > 

In all _E-integrals the numerical integration is done up to E = i5 max and the asymptotic con- 
tribution from i?max up to oo is evaluated analytically and added to the numerical result. 
Similarly, an upper limit for the a-integration is determined by requiring that the it-integral 
over the ubiquitous exponential e(tt, a) in Eq. 14.49H 

l-« "\ 1 2 2, s 1 (M/2fa 2 

~ 2 m ^ (<J) = 2^>) 



K(x, y) = I du exp ( -x — yu) , x = -m n (a) , y = - — 7T r^ — (D.3) 



should be smaller than e = 3.1 • 10 7 . Using the large-cr limit of this expression from Eqs. 
(82) - (84) in ref. |26l we therefore have to take 

(D.4) 



~~ rnM 

Also, in all <r-integrals the numerical integration is done up to a — <r m ax and the asymptotic 
contribution from <j max till oo is evaluated analytically and added to the numerical result. 
This strategy works very well as we know the large (E, a)-behaviour of profile functions and 
pseudotimes (see appendix C). 
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The number of gaussian points n g (typically 72) and subdivisions n e (typically 3 — 4 ) 
should, of course, be adapted to the values of -Emax, fma X : ideally at least one point should be 
on each oscillation of the trigonometric function sin 2 (Ea /2) = ( 1 — cos (Ea) ) / 2 for all values of 
E,a but experience showed that £ m axffmax/(Tfi 9 n e ) ~ 3 — 4 was sufficient. 

In principle, the it-integral for K(x, y) in Eq. 1D.3I can be expressed in terms of a (non- 
standard) special function, i.e. either by a special Shkarofsky plasma dispersion function (Eq. 
(26) in ref. |73| 1 or by an incomplete MacDonald function of first order in Bessel form (as defined 
in Eq. (III. 1.8) of ref. |74l~). However, by doing so nothing is gained for practical computations. 
Therefore, for simplicity and ease of implementation, we have used a direct gaussian integration 
of the defining integral ljD.30 with the same number of gaussian points and subdivisions as for 
the solution of the variational equations. Of course, for m — this is not necessary as the 
u-integration can be done analytically. 

The problem of numerical computation of integrals with oscillating integrands is also 
present when one tries to evaluate the virial theorem in the form of Eq. 15.131 . Here the 
derivatives of the pseudotimes 



da 



dE 



+ (-) 
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A-(E) v ' A+(E) 



(D.5) 

E 



are needed which have oscillating integrands with less damping at large E than the pseudotimes 
themselves. This may lead to a loss of accuracy in the numerical evaluation. It can be shown 
that for theories without derivative interactions - like the WC model but not for QED - the 
derivative term in the virial theorem 15.131 may be eliminated completely. However, Tabled 
and show that below the critical coupling such a form of the virial theorem is not needed as 
the virial check is satisfied within a factor of two of the relative deviation which measures how 
well the variational equations are fulfilled. 

If one solves the coupled variational equations by iteration there is a crucial difference now com- 
pared to the one-body case in that one cannot start with A+(E) = 1 as this would correspond 
to an unbound system. Instead we have used the nonrelativistic solution A+(E) = 1 +u> 2 /E 2 
with the Coulombic value I4.6UH for d = 4. Obviously, the Yukawa value would be a better 
starting point although that doesn't really matter since one has to iterate many times anyway. 
Of course, the singular behaviour of the profile function A+(E) at E = has to be treated with 
care; this is done by updating the coefficient u vaT from Eq. 1C.3I each time during the iteration. 
To shorten the iteration one can use the A-value from the one-body case or the previous solu- 
tion as starting point. As a measure of convergence during iteration we have used the relative 
deviation 



Ai(a;) 



x 



(i) 



(D.6) 



for each variational parameter x = A, tj var and the maximal relative deviation of the variational 
functions x = m&x n [A±(E n ) , /U 2 ^ ; 12 (<x„)] at each gaussian point. From these relative deviations 
a global measure of convergence 

A, = max[Ai(A),A ),A l (A±),A l (/i 2 1/12 )] (D.7) 

was constructed and monitored during iterations. The updating after each iteration was per- 
formed with a prescription given by Hauck et al. {75 1 which suppresses oscillating values from 
one iteration to the next by using an exponential weigthing of old and new values: 

_ e -AiM) x (i+l) + x (i) a .(i+l) < (D . 8) 

This is most effective in the early stages of the iteration when the relative deviation Ai(x) is 
still large. 

The iteration was stopped when the global measure Ai was less than a prescribed value 
for which typically 10~ 3 was taken. Although this is far away from the 10~ 5 - accuracy which 
was easily obtained in the one-body case, it should be sufficient for the present purposes. 
Given the much more demanding numerical problems of the two-body case we believe that any 
improvement would need much more computing power (i.e. more integration points) and/or 
more efficient algorithms (761 177| to solve the nonlinear integral equations. 
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Appendix E: Critical oupling Constant and Width: Analytic Approx- 
imation 

Here we show how one obtains the critical coupling from the simple ansatz A- (E) = 1, A+(E) = 
1 + lj' 2 /E 2 , keeping A and w as variational parameters. For this choice and ke = 1 one easily 
finds the kinetic terms and pseudotimes as 



fl_ = , fi+ = to , /Ltij(er) 

Furthermore we assume weak binding 



a 

2 + 2lo 



l - [l + (-y e -^] (7 = 1,2). (E.l) 



(9/2) 2 



M 2 



<C 1 



(E.2) 



and that cua <C 1 in all integrals over the proper time. Then we may take approximately 



and obtain after performing the <r-integration 



V n ~ — M 2 
2n 



du In 



1 + 



- - 



(\q/2) 2 u 2 



1 - u 



via 



-.M 2 



2tt 



du 



2^l 



exp 



m 1 — u 
2u> u 



(E.3) 

(E.4) 
(E.5) 



Equation (IE.4I is identical with Eq. (46) in ref. [II] if we use the weak binding approximation 
q/2 = M. Mano's equation then approximately reads 



M = (2A-A 



(2A-A 2 ) (|) _ w _2(Vii+Via) 



(E.6) 



We now vary w.r.t. A and to and obtain (for the case m = where the it- integrals can be 
performed analytically) 



I - Q M 2 

A ~ + K(\q/2) 2 

(Za) 2 M 4 
W 2tt (Ag/2) 2 



l + Z 



2(Xq/2) 2 



(E.7) 
(E.8) 



Note that Eq. ljE.80 reduces to the d = 4-result in Eq. I)4.6()|l if we replace Ag/2 by M. Here, 
however we keep the A-dependence which in the one-body case was essential for probing the 
instability of the WC model. Inserting Eq. <E.81 into the A-Eq. <IE.7t and taking q/2 ~ M now 
gives a quartic equation for the variational parameter 



A 4 -A 3 +°A 2 + ^ = 
7r 2n 



(E.9) 



which for Z = reduces to the approximate one-body equation (149) in ref. [I]. To find the 
critical value of the coupling constant where this equation ceases to have real solutions we do 
not have to use the cumbersome explicit solutions; rather we observe that branching into the 
complex plane occurs when d\/da = 00. Therefore we may differentiate Eq. 1E.9I w.r.t. to A 
and setting da/dX = we obtain 



4A 2 rit - 3A crit + 2^=0. 



(E.10) 



If this is put back into Eq. <E.9t we can solve for 



Acrit = t ^— - — \l-3z + Vl + 3z\ 



M-z 



where 



2-kZ 2 



(E.ll) 
(E.12) 
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and obtain the critical coupling as 

* 1-9* + (1 + 3*)"* (E . 13) 



Oicrit = ~ 
O 



{1-zf 



By simple algebraic manipulations this can be brought into the form lb. 181 . 

To obtain an analytic estimate for the width of the ground state above the critical coupling 
one allows the product C, = \q/2 =: £o exp(— i\) to become complex. The quartic equation 
HE. 91 then reads 

(E.14) 



a M (ZaY M 
H h ^ - 

7T C 2 27T C 4 



If this and the variational solution are inserted into the approximation ljE.60 to Mano's 

equation we obtain 



,2 ,-,2 2o< 2 3(Za) M Q 2 
C = M M ^ — —7T- + —M 



TV 



7T 



du In 



m 2 1 — u 



(E.15) 



which is the generalization of Eq. (48) in ref. [II]. For the imaginary part of this equation it is 
possible to set immediately the meson mass to zero 18 with the result 

3(Za) 2 M 4 



Co sin 2x 



— M 2 X 
n 



(E.16) 



(compare with Eq. (49) in ref. [II]). Equation IE. 1611 allows the modulus £o to be determined 
in terms of the phase x- 



r 2 



—M 



[l+ S 2 ( X )] 



+ (E.17) 



7T 2S(x) 

For Z = this reduces to Eq. (49) in ref. [II]. The width is obtained from Eqs. 16.231 and HE. 141 



T (Za) 2 M 4 a M 2 



i + 



* Co 2 



(Za) 2 M 4 
H — cos3x 



by taking real and imaginary parts: 
M = ( cosx 
T = 4Cosin X 

With the help of Eq. IjE.lBfl the last equation can be brought into the form 

r 



_ 1 — 1 



27T Q 

2(Za) 2 M 4 
Co 



sin 3 X 



(E.18) 

(E.19) 
(E.20) 

(E.21) 

(E.22) 
(E.23) 

7T Z 

is the value of the modulus at the critical coupling obtained from Eq. ljE,170 by setting \ — 0. 
The parameter z has been defined in Eq. IE. 121 . 

It remains to determine x(°0- As in ref. (II) one obtains the inverse relation by combining 
Eq. l|E~T9t with Eq. ijETfjl . After some algebra this gives 



4Cosin X [l_5(x)] + ^^^ [3sin X S(x)-sin3x] 
77 so 



which shows that close to the critical coupling (x — + 0) the width behaves as 



Here 



t 2 - 

Sent 



8 ^ ■ 

^ Sent 



Qcrit 



1 + 



3(Za clit ) 2 M 4 
^ C c 4 rit 



M - [l+VT + 3i] 



a = 2^(x)Kz,x)f 

• {cos X [l + S(x)] [w(z, x )} 2 + zS 2 ( X ) [cos3x-3 5(x)cos X ]}" 



(E.24) 



18 This can be also done for the real part if we eliminate M from Mano's one-body equation, e.g. 
by setting Z = in Eq. HE. 151 . 
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where w(z, \) = 1 + v/l + 3zS(x). This result is the generalization of Eq. (53) in ref. [II] to 
which it exactly reduces for Z = z — 0. Expanding a around the critical value in powers of \ 
yields 

a = ttcrit [l + aax 2 + ---] (E.25) 

where 

a 2 = 1 + 2 ^ + ^ T ^ . (E.26) 

1 + z + VI + 3z 

Finally instead of Eq. (57) in ref. [II] one gets from Eq. (IE. 2211 the following width for the 
ground state 



T ~ -2M 

3 



a — Qcrit 
«crit 



3/2 



/corr {Z) + . 



(E.27) 



where the explicit form of the correction factor is given in Eq. 16.26H . 
Appendix F: Weak-Binding Limit 

Here we derive the weak-coupling limit for the binding energy in the case of massless pions 
(which is the original Wick-Cutkosky model). Setting m = allows us to perform the calcula- 
tions analytically. First this is done for the variational approximation and then in the potential 
version of effective field theory appropriate for the present model 78 . 

F.l Variational Calculation 



For the variational approach a convenient starting point is the Feynman-Hellmann theorem 
115. 1711 for which one only has to evaluate 



V 12 = -« M2 



tt (Ag/2)2 



da — 



exp 



(A?/2)V 



2^ 12 {a) 



(F.l) 



Again, we have choosen to work in the reparametrization "gauge" Kb = 1- From Eq. ljE.10 
we infer that approximately u 2 2 (a) ~ 1/u) + uia 2 /A + . . . with u) — O {(Za) 2 ) (see Eq. {Ej 
Therefore we insert the expansion 



/"12(f) = uo + u 2 a + . 
into Eq. I|F,10 . A simple calculation then gives 



, U2 < 110 



(F.2) 



V12 = 



Za M 



2tt A<?/2 y/uo~ 



1-i- 



11 2 



2 (A 9 /2) : 



+ . 



(F.3) 



We have to determine the coefficients uo,u 2 from the variational equations. This we can do 
perturbatively except that the low-E behaviour A+(E) — ► uj^/E 2 needs to be included to all 
orders to generate a bound state. Thus from Eq. 15. 121 and Eq. 15.51 1 we have 



1/0 



2 

7T 
1 

^var 

1 

^var 







dE { E 2 + u&r + 2 / da [5Vu - 8V12] sin 2 (Ea/2) 
4 



da [SV11 - 8V12} 



dE si » 2(J ^!l + . 



(u>^ T + E 2 ) 2 



da a 2 {SVu - SVn) 



(F.4) 



Here we have used u> vaT a <C 1 and the abbreviation SV\j 
is inserted into the equation defining w^r 

Za M 2 1 



SVxj/SnlAa). If the result JFtJ 



daSVw = 2 



du 



we obtain 



,V2 



Za M 



2^ Xq/2 



— [ da a 2 
4 Jo 



2^ Ag/2 u 3/2 



1-i 



"2 



2 (Ag/2)2 



-3/2 



1-i 



1/2 



2 (Ag/2)2 



(F.5) 



(F.6) 
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The correction terms may be safely expanded and the approximations A ~ 1, q/2 ~ M, /Mi! 
a, ftia — l/ w > u '2 — w /4, w — (Za) 2 /(27r) used without impunity. Thus we obtain 



(Za) 2 



M 



1 + ^ 



5 (Za) 



2tt (Ag/2) 2 

Note that the RHS of Eq. l|F.5jl is just — V12/U0. Hence 

{Za) 2 M 4 



4tt 16 7T 



+ . 



Via 



2 

-WO Wvar 



2tt (Ag/2) 2 



1+ 2^ 



4% 



+ 



(F.7) 



(F.8) 



It remains to determine the variational parameter A up to order a, (Za) 2 . 

This can be done by using the variational Eq. 15.31 with the zeroth-order approximations 
for the pseudotimes or simply by expanding Eq. IIE.71 which is correct to that order. In both 
cases one obtains 

a | (Za) 2 ' 

7T 



A = 1 - 



2tt 



+ • 



(F.9) 



which now allows application of the Feynman-Hellmann theorem 15.171 



q\ 2 d 



dZ \2 



Za , ,4 
M 

TT 



„ 7 a 5 Z 2 a 2 



2tt 4 1 2 



Integrating on both sides and using q/2 = M + eo gives 



(M + eo) 4 = const. — 



(Za) 2 



M 



1+ 2n 



8 IT 2 



(F.10) 



(F.ll) 



The integration constant must equal M 4 because the two particles are unbound for Z = 0. 
Therefore the ground-state binding energy has the weak-coupling expansion 



M 

60 = 2 



(Za) 



7 a . 

1+ 2n ] - 



(Za) 4 



+ ■ 



(F.12) 



F.2 Effective Field-Theory Calculation 

For simplicity, we consider here not the Lagrangian 11.21 but a system of identical nucleons 
described by the Lagrangian 



C = & (-d 2 -M 2 +g' X )$ + £o(x) 



(F.13) 



where the last term is the free Lagrangian for the mesons. In the nonrelativistic limit ( M — > 00) 
we make the ansatz 



$(x) = 



1 



I2M 



-iMx 



(x, Xo = t) 



(F.14) 



to describe particles (anti-particles would have a different sign in the phase factor and are 
omitted as explicit degrees of freedom in the following). This leads to 



C 



A g' 
ldt+ 2M + 2M X ' 



1 d 2 I 4> + c ( x ) 



2M 



(F.15) 



The action is the 4-dimensional integral over the Lagrangian (density) so that an integration 
by parts brings the last term in the brackets into the form (fr(f)/(2M) . Using the equation of 
motion for the nonrelativistic field 6 



one obtains 



Therefore 



(~2iMd t -A + d 2 )<f>(x,t) = fl 'x0(x,i) 



(*.*) - 7^ (A + fl'x) <Kx,t) 



c. 



2M 



(2Mf 



6+ (A + g'x) 



(F.16) 
(F.17) 
(F.18) 
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is the l/M 3 -correction to the leading nonrelativistic Lagrangian. Its form is in agreement with 
Eq. (3.6) or Eq. (4.12) in ref. [79 and can also be easily obtained from the Foldy-Wouthusen- 
Tani Hamiltonian <j>* \J M 2 — A — g'\ </> by an expansion in inverse powers of the heavy mass 
M. Note that a "seagull" term 4>* r x' 2 4> appears in the nonrelativistic field theory. 

Since one has changed the high-energy behaviour of the theory the correct effective La- 
grangian (for particles) is 

+ ftp ^ i 9 ' 02 P X + + 9 ' 2 ° 3x2 + 9 ' dl (Ax) ] + ■ ■ ■ + £o(x) ■ (F - 19) 
This should be accurate up to order Mv 4 in the (velocity) counting rules 19 in which 
dt = O (M« 2 ) (nonrelativistic energy) 

A — O (M 2 ?; 2 ) = O (M« 2 ) (nonrelativistic kinetic energy) 

-^-^ = V A-Koi = O (y l ^J (velocity in Bohr orbit = a) 

X = O (Mv z/2 ) =► ^ X = O {Mv 2 ) (potential energy) . (F.20) 

In Eq. IF. 1911 the "Wilson-coefficients" a = 1 + 0(g' 2 ) encode the missing high-energy (short- 
distance) information. They are obtained by "matching", i.e. by comparing physical amplitudes 
in the full and the effective theory. The coefficient di is expected to be 0{g' 2 ) since it is not 
present in the tree-level calculations and this is confirmed by the explicit calculation below. 
Note that the kinetic terms A/(2M ), A 2 /(8M 3 ) are not renormalized since they reflect the 
exact Lorentz symmetry of the underlying relativistic theory. 

We determine the Wilson coefficients ci,C2,di by evaluating the meson- nucleon scattering 
amplitude in both theories including the lowest-order radiative corrections. To do that we also 
need the residue Z r of the nucleon 2-point function 



p 2 - Ml + E(p 2 ) p 2 - M 2 

where M 2 = M 2 + E(Af 2 ) and 

9S(p 2 ) 



Z" 1 = 1 + 



dp 2 



(F.21) 



(F.22) 



J 2 = A/ 2 



Of course, in lowest order we have Zi°^ = 1 . Using Muta's conventions |80| we obtain for the 



1-loop self-energy (depicted in fig. 



v(l)f 2n /2 4-d 

^ (p ) = g v 



d d k 1 1 

(2iv) d i M 2 -(p-k) 2 m 2 - k 2 



= , 2 V{2-d/2) r 1 

9 (47r) d / 2 J [M 2 x + m 2 (l-x)-p 2 x{l-x)] 2 - d / 2 ' { ' 

Here we have introduced standard Feynman parameters and performed the diverging momen- 
tum integral in d dimensions. In this order we may replace the bare mass Mo by the physical 
mass M, differentiate w.r.t. p 2 , expand and find 

7 (D _ „/2, 4-d r(3 - d/2) f 1 ^ x(l-x) 



= -9 V / dx 



2 



(4yr) d / 2 J [M 2 x 2 + m 2 (l- x)f~ d / 
d ~ ^ fl (F.24) 



16tt 2 J M 2 x 2 + > 

19 Since v = v/c the counting is equivalent to the expansion in appendix lA.ll but the procedure 
is a systematic one which also allows to include loop effects. 
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which is UV-finite but IR-divergent for m = 0. 

The tree-level truncated nucleon-meson (2,l)-point function (multiplied by —i) is just the 
coupling constant g' . The one- loop correction to this vertex function is shown in fig. [|)p and 
leads to 



d 4 k 



1 



1 



(2tt)H M 2 - (p - k) 2 A/ 2 - (p - k + q) 2 m 2 - k 2 
For g^Owe obtain by standard techniques 



(F.25) 



T W (p,q^0) = 



16tt 2 



dx 



M 2 x + m 2 (l — x) ~ p 2 x(l — x] 



+ 0(g' 3 q 2 ) 



(F.26) 



The physical amplitude for meson-nucleon scattering has two external nucleon legs and therefore 
requires twice a wavefunction renormalization constant \fZ^- applied to the truncated (2,l)-point 
function. Hence 



T(q^0) = 



9cH 



Z r T(p,q^0) 

12 
l07T 



= • ffeff 



p 2 =M 2 
1 

dx 

o 



+ 0(q 4 



M 2 x 2 +m 2 (l - x) 



(F.27) 
(F.28) 



Here (r 2 ) is the root-mean square radius of the nucleon due to radiative corrections. 

The meson-nucleon scattering amplitude in the nonrelativistic theory described by the 
Lagrangian l|F~19)l reads 



1 (P ' q) = < ! Cl 2M + (2MF 



C 2 • 



-p 2 - (p + q) 2 ) - <W] 



\J2E V 2E P 



(F.29) 



where the square-root factor is due to the different normalization of relativistic and nonrelativis- 
tic single-particle states 
we obtain 



Expanding E p = \J M 2 + p 2 and E p +q for low three-momenta 



Tnonrcl. / 
(P 



0, q — > 0) = cig' + (ci - c 2 ) 



, p 2 + (p + q) 2 

7 4M 2 



d.2 g 



4M 2 



+ . 



(F.30) 



Comparing with Eq. 1F.27I we find for massless mesons 



ci = c 2 = — - = H 



(F.31) 



and (q 2 



-q 2 when all momenta are small) 



dl = ^M!»f(r 2 ) = 0(g' 2 ) 
3 ff' 



(F.32) 



The coefficients ci , C2 account for an enhancement of the effective meson-nucleon coupling 
constant through loop effects similar as in Schwinger's famous determination of the anomalous 
magnetic moment of the electron in QED. The coefficient d\ describes the finite extension of 
nucleons due to the meson cloud and therefore is infrared divergent. This can also be seen from 
its explicit lowest-order perturbative expression as given in Eq. (60) of ref. |2(i| but we do not 
need it in the following. 

In the next step one has to integrate out the mesons in order to reduce the relativistic 
bound-state problem to a quantum-mechanical one. This we do by solving the equation of 
motion for the meson field (operator) obtained from varying Eq. !F.I9t 



(d + m 2 ) xo 



2M 



(F.33) 



Since in Eq. 1F.3H we have determined Ci only to one-loop order, i.e. to order v in the counting 
rules, we have retained only the leading term and for consistency have neglected the 0(v 2 ) 
corrections in the effective Lagrangian. Furthermore, we replace g — > Zg to indicate that the 
mesons are from the other particle (self-energy effects vanish in the effective nonrelativistic 
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theory if dimensional regularisation is used 
present the solution of Eq. ljF.330 is 



Zg' 



For the case that no external mesons are 



Zg' 



d 2 + : 



— — -c\ Id y ( x 
2M 1 y X 



1 



-A + m 



y) <j>Hy)<P(y) 



(F.34) 



According to the counting rules, the time derivative in the d'Alembertian is also suppressed 
by a factor 0(v 2 ) compared to the Laplacian. Equation IF. 3411 is therefore correct up to order 
Mv 5/I2 and after substitution into Eq. (IF. 1911 gives the following effective (potential) Lagrangian 

£ c p ° 4 = ^(x.i) ji0 t + A_I J d 3 V cj ) \y,t)V eS (^-y)4>(y,t) + 0(Mv i )^ #x,t) (F.35) 



where 



Kff (x - y) = -Cl 



(z£ 

\2M 



1 



-A + m 2 



2 Za _ m | X _y| 

y ) = _c i 1 r e 



(F.36) 



is the usual Yukawa potential enhanced by a factor c\ from radiative corrections. In the Coulom- 
bic case (m = 0) this means that the leading term for the binding energy is 



M (Za) 2 
~2 2(n + 1} 



4 
Cl- 



M (Za) 2 
~2 2(n + 1)* 



[l + 4^+0((Za) 2 a 2 ,(Za) 4 



,n = 0,1,2. 



(F.37) 
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